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INTE    .'TOW 

An  important  operating  mode  of  deep  submersibles  is  the  ascending 
and  descending  mode.  Unlike  large  high  speed  military  submarines  V7hich 
ascend  and  descend  dynamically  by  use  of  a  combination  of  speed  and  append- 
age deflection,  deep  submersibles  depend  partially  and  sometimes  entirely, 
as  in  the  case  of  the  TRIESTE,  upon  buoyant  ascent  and  descent.  Of  the  as- 
cending and  descending  operations,  the  most  critical  is  during  an  emergency 
ascent  in  which  most  or  all  of  the  jettisonable  weights  are  removed  from  the 
vehicle  in  order  to  obtain  as  rapid  an  ascent  as  possible. 

The  majority  of  deep  submersibles  fall  into  a  class  of  vehicles 
characterized  by  low  speed  and  minimal  appendage  streamlining ■  They  are  de~ 
signed  for  underwater  research  and  rescue  missions  which  do  not,  except  for 
the  travel  mode  of  the  Deep  Submergence  Rescue  Vehicle,  require  large 
speeds.  They  are  equipped  with  many  and  diverse  appendages  for  accomplish" 
ing  their  missions  such  as  manipulator  arms,  lights,  propulsion  motors, 
sampling  equipment,  TV  cameras  and  mating  skirts,  all  of  which  effect,  to 
varying  degrees,  the  streamlining  and  hydrodynamics' of  these  vehicles. 

In  order  to  study  the  motions  and  ascent  trajectories  of  sub- 
mersibles in  normal  or  emergency  ascent  without  model  or  full  scale  tests, 
it  is  necessary  to  develop  a  mathematical  model  of  the  vehicle  which  can 
accomodate  the  effects  of  varying  mass  and  center  of  mass. 

This  work  is  an  attempt  to  develop  the  mathematical  model  using 
a  set  of  second  order  equations  of  motion,  and  solve  these  equations  by 
use  of  a  high  speed  digital  computer. 

The  primary  sources  of  material  for  this  work  are  the  works  of 
Abkowitz,  Dogan,  Gertl  ,3   b  and  Strumpf,  but  the  secondary 


sources,  acknowledged  or  not,  are  no  less  important  ot  the  development  of 
this  wo:  . 

The  problem  of  determining  the  motions  of  a  totally  submerged 
body  has  been  treated  by  many  authors,  however,  there  appear  to  be  only 
two  basic  means  for  arriving  at  these  equations.  These  are  by  energy 
methods  and  by  vector  calculus. 

The  first  method  is  that  used  by  Sir  Horace  Lamb  in  his  book 
"Hydrodynamics"  which  first  appeared  in  1879  (see  ref  1).  The  enorgy 
method  is  based  upon  the  existence  of  a  single  valued  velocity-potential 
which  implies  that  tho  motion  of  the  fluid  is  entirely  due  to  that  of  the 
submersible,  and  is  therefore  irrotational  and  acyclic.  This  method  leads 
to  a  completely  general  set  of  equations  representing  the  rigid  body  dynam- 
ics, but  the  hydrodynamic  forces  that  result  represent  only  the  so  call',  i 
added  mass  or  added  inertia  effects  and  not  all  the  forces  acting  on  the 
body.  Of  the  forces  that  are  missing  the  most  important  are  those  due  to 
circulation,  separation  and  vortex  shedding, 

Tho  second  method  of  developing  the  equations  of  motion  is  to 
combine  a  vector  expansion  of  Newton's  laws  of  motion,  expressed  as  follows: 

F  =  £j:'rr  (Momentum)  -   tho  vector  force  acting  on  the  body        (l) 

M  =  -rr  (Angular  momentum)  =  the  vector  moment  acting  on  the 

body  (2) 

with  a  Taylor  series  expansion  of  the  hydrodynamic  forces  and  moments  ex- 
pressed as: 

foments  ~  f  (Properties  of  body;  Properties  of  motion, 

Properties  of  ,  (3) 


which  reduces  to 


F 


—  =  f  (Properties  of  notion)  (4) 

for  a  particular  vehicle  in  a  particular  fluid.  This  method  has  been  used 
by  A.  Strurapf  in  i960  in  developing  his  "Equations  of  liotion  of  a  Sub- 
merged Body  with  Varying  Mass"  (see  ref  2)  and  by  Prof,  M,  A,  Abkowitz  in 
19^9  in  his  lecture  course  and  notes  on  "The  Dynamical  Stability  of  Sub- 
marines" (see  ref's  3  and  10),  It  too  results  in  a  completely  general  set 
of  equations  describing  the  rigid  body  dynamics,  but  its  greatest  asset 
lies  in  the  generality  of  the  hydrodynamic.  forces  and  moments.  The  hydro- 
dynamic  forces  and  moments  obtained  from  a  Taylor  series  expansion  in- 
clude not  only  all  the  added  mass  effects  but  also  the  circulation,  sepa- 
ration and  vortex  shedding  effects.  The  only  limitation  on  this  metfr 
is  the  availability  of  theoretical  or  experimental  data  to  use  in  the  equa- 
tions. 

From  the  preceeding  discussion  and  equations  (i  through  *0,  it  is 
obvious  that  the  derivation  of  the  rigid  body  motion  may  be  completely  sep- 
arated from  the  problem  of  developing  a  suitable  form  for  the  hydrodynamic 
effects.  Therefore,  the  two  methods  of  development  may  be  broken  down  and 
part  of  each  used. 

This  work  shall  attempt  to  take  advantage  of  the  simplicity  of 
the  energy  development  of  the  rigid  body  dynamics,  while  retaining  the  gen- 
erality provided  by  the  Taylor  series  expansion  of  the  hydrodynamic  forces 
and  moments. 

The  solution  of  the  equations  of  motion  that  will  be  developed 
employs  a  technique  suggested  in  196*1  by  C.  Parissis  in  his  solution  c f 
linearised  equations  of  motion  for  heave  and  pitch  of  a  ship  (see  ref  H) 


and  again  in  19&7  by  *-■•  M,  McCloskey  for  tho  digital  simu]     i  of  the  D.SRV 
control  system  and  autopilot  (see  ref  5)<> 

Tho  coefficients  used  to  test  the  solution  of  the  equations  of 
notion  come  from  a  series  of  model  experiments  conducted  on  tho  DSHV,  The 
results  of  the  computer  solution  will  not,  however,  be  the  ascent  trajectory 
for  the  B5RV  since  there  is  an  assumption,  in  this  work,  of  zero  propulsion 
forces  in  the  development  of  the  equations  of  motion. 

The  assumptions  and  the  statement  of  the  problem  are  given  in  Chap- 
ter 1,  the  equations  of  motion  are  derived  in  Chapter  II  and  a  method  for 
solving  the  equations  of  motion  is  developed  in  Chapter  III. 


CHAPTER  J 
FORMULATION  OF  THE  PROBLEM 

The  problem  involves  a  deep  submersible  vehicle  of  given  geometry 
and  dimensions  which  is  in  free  ascent  through  a  stationary  fluid. 

I—  1  Assumptions 

The  vehicle  is  assumed  to  be  a  rigid  body  with  no  elastic  defor- 
mations of  a  vibratory  nature*  It  has  six  degrees  of  freedom,  three  trans- 
lational  and  three  rotational.  We  shall  bo  interested  in  motions  in  both 
the  horizontal  and  vortical  planes  of  motion.  Velocities  arc  small •  Hydro- 
dynamic  effects  of  second  order  in  acceleration  shall  bo  considered  negligi- 
ble. The  only  body  symmetry  is  port  and  starboard. 

I  -  1.1 

The  neglecting  of  vibratory  motions  is  reasonable  since  the  fre- 
quencies of  vibration  of  the  hull  acting  as  an  elastic  body  are  of  different 
orders  of  magnitude  than  the  frequencies  of  motion  and  do  not  excite  the 
latter.  Elastic  deformations  due  to  vehicle  compress ability  must  be  in- 
cluded since  they  directly  affect  the  buoyancy  of  the  vehicle* 

I  -  1.2 

The  interest  in  both  pirn  is   of  motion  is  duo  to  the  desire  for 
as  general  a  sot  of  equations  as  possible. 

I  -  1.3 

The  assumption  of  small  velocities  is  realistic  ."  {'  e  as- 

cent velocities  and  open  ocean  currents  are  in  general  of  order  less  ' 
ten  knots. 


"I  »  1.4 

The  usg  of  a  sot  o.f  second  order  equations  of  motion  i*  d  mod 
necessary  to  appropriately  describe  the  hydrodynamic  cross-coupling  that 
takes  place  in  a  problem  such  as  this.  The  second  order  acceleration  ef- 
fects are  assumed  to  be  zero  on  the  bas.".r>  of  potential  theory  (sec  rof  l). 
This,  however,  has  not  been  experimentally  verified 0 

1-1.5 

Neglecting  as symmetries  due  to  relatively  small  appendages  which 
are  not  control  surfaces,  there  are  few  vehicles  which  operate  in  the  ocean 
environment  which  do  not  possess  port  and  starboard  symmetry.  For  this  rea- 
son the  assumption  of  port  and  starboard,  symmetry  has  been  made. 

1-2  Initial  and  Other  C  one! itions 

The  vehicle  is  to  be  initially  at  rest  relative  to  the  inertial, 
earth  fixed,  axes.  The  control  effectors,  propellers,  thrusters,  rudders, 
dive  planes,  etc.,  are  inoperable  and/or  in  neutral  position.  The  driving 
force  for  the  vehicle  shall  be  a  decrease  in  weight  due  to  jettisoning  of 
ballast  or  an  increase  in  volume  due  to  vehicle  decompression  as  it  rises. 

These  conditions,  though  arbitrary,  serve  to  restrict  the  problem 
to  one  of  manageable  size. 


CHAPTER  II 
EQUATIONS  OF  MOTION 

II  -  1  General^ 

The  derivation  of  the  equations  of  motion  for  a  rigid  body  in  six 
degrees  of  freedom  with  a  varying  mass  and  center  of  mass  follows  closely 
that  of  Dr.  Pierre  Dog an  in  reference  (5)1  which  is  based  primarily  upon 
the  development  of  the  equations  of  motion  by  Lamb  in  reference  (l).  The 
derivation  of  the  hydrodyhamic  force  equations  follows  that  of  Professor 
Martin  Abkowitz  in  reference  (3)  except  that,  in  this  work,  the  second 
order  terms  are  retained.  The  forces  due  to  gravity,  are  determined  using 
the  transformations  set  forth  in  reference  (6)  and  the  angular  velocity 
transformations  given  by  Professor  Abkowitz  in  reference  (7). 

The  axis  systems  necessary  to  describe  the  motions  of  a  body  and 
its  trajectory  through  a  fluid  include  a  body  fixed  system  and  an  inertial 
system.  » 

H  -  2  Axis  Systems 

The  equations  of  motion  for  an  ascending  vehicle  must  be  written 
in  an  earth  fixed  axis  system  in  order  to  determine  the  trajectory  of  the 
vehicle  relative  to  some  fixed  point.  An  additional  body  fixed  system  is 
required  in  order  to  describe  the  hydrodynamic  interactions  between  the 
vehicle  and  the  water. 

II  -  2.1  Inortial  Axis  System 

The  right-handed  earth  fixed  axis  system,  A.,  (see  figure  11,1 ), 
is  assumed  to  be  an  inertial  syst  n  for  the  i   son  that  the  accelerations 


AG      BOD/AXES 


/\  EARTH     FIXED 

AXES 


Figure  1  Coord 
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of  a  point  on  the  surface  of  the  earth  are  an  order  of  i  ,     le  sraal] 
than  those  which  are  of  importance  to  tho  motions  of  the  vehicle.  A_  is  an 
orthogonal  set  of  axes  fixed  relative  to  the  earth  such  that  components  x^ 
and  y„   are  in  a  horizontal  plane,  and  the  z,-,  axis  is  vertical  and  directed 

ti  £j 

downwards • 

II  -  2.2  Body  Axis  Sys-' 

The  right-handed  body  axis  system,  A_  (see  figure  II.l),  is  fixed 
in  tho  vehicle  such  that  advantage  is  taken  of  the  assumed  principal  plane 
of  symmetry  by  placing  tho  origin  of  the  system  in  this  plane.  The  axes  of 
this  system  are: 

x  axis  -  the  longitudinal  axis,  directed  from  the  after  to  the 

forward  end  of  the  vehicle, 
y  axis  ~  the  transverse  axis,  directed  to  starboaid, 
z   axis  -  tho  normal  axis,  directed  from  top  to  bottom  (deck  to 

keel) . 

The  xz  plane  is  the  assumed  principal  plane  of  symmetry. 

11  -  2.3  Body  Axes  Orientation 

Angular  displacements  of  Ag  relative  to  A.-,  are  specified  by  a  set 
of  modified  'Euler  angles'  which  are  taken  as  positive  in  the  sense  of 
rotation  of  a  right-handed  screw  advancing  in  the  positive  direction  of  the 
axis  of  rotation. 

The  orientation  of  Ap  relative  to  /  is  described  in  terms  of  a 
roll  angle  $5,  a  pitch  angle  G  and  a  yaw  angle  "f.  Before  defining  these 
angles,  an  order  of  rotations  must  be  chosen  since  finite  rotations  are  not 
true  vector  quantities  and  do  not  obey  the  rules  for  adding  vectors.  T] 
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order  chosen  here  conforms  to  that  set  forth  in  reference  (6)  which  is: 

(1)  rotate  about  the  initial  z   =  zp  axis  through  an  angle  of 
yaw  f, 

(2)  rotate  about  the  new  position  of  the  y  -  y,  axis  through  an 
angle  of  pitch  G, 

(3)  finally  rotate  about  the  new  position  of  the  x  -  x  axis 
through  an  angle  of  roll  £<, 

In  accordance  irith.  the  order  of  rotations  above  wo  have  the 
f ollowing  definitions : 

G  -  the  angle  of  pitch j  the  angle  of  elevation  of  the  x  axis}  the 
angle  between  the  x  axis  and  the  horizontal  plane  x^y^, 

"jr-   the  angle  of  yaw?  the  angle  from  the  vertical  plane  x-,z„  to 
the  vertical  plane  xz^, 

<}>  -   the  angle  of  rollj  the  angle  from  the  vertical  xzr  plane  to 
the  principal  plane  of  symmetry  xze 


The  successive  rotations  required  to  specify  the  orientation  oj 
the  body  axes  relative  to  the  earth  fixed  axes  can  be  described  by  three 
orthogonal  matrices  [yl  ,  (_CJ   ,  fjM  : 


cos  f 

-sin  ^r 

0 

cos  0 


sin  if 
cos  y 
0 


0 


0 
0 

1 


■sin  e 


0     1       G 
sin  C     0      cos  C 


a = 


0 

CO 

0 

-sin  <f> 

;in  o* 


COG  J; 


Unit  vectors  in  AR  and  A,,  are  related  by  an  orthogonal  matrix: 


*B 


E 


- 

- 

V 

'h 

Jb 

=s  T 

k 

kB 

X 

where  T^  is  the  product  of  the  three  orthogonal  matrices  defining  the 
rotations : 


cos  6  cos  f 

•cos  <f>  sin^ 
+sin  $  sin  G  cos  V" 

sin  ^  sin  ^ 
•fsin  ©  cos  §  cosf- 


cos  0  sin  ir 

cos  fJ  cos^ 

+sin  ^  sin  0  sin  "fr 

-sin  ^i  cos  y* 

+cos  $  sin  6  sinV'" 


-sin  G 
sin  p*  cos  G 

cos  <f>   cos  C 


Velocities  in  AR  and  A^  are  also  related  by  the  transformation 


matrix  T^: 


-    -1 

u 

- 

• 

V 

=  TB 

^e 
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The  vehicle  angular  velocities  in  Ag  and  A„  are,  however,  related 
by  a  non- orthogonal  matrix,  which  is  the  sum  of  three  components  along  t! 
ZE'  yn  an<^  x  axos  0:f  mao^itude  -yv.  0  and  $. 
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q 


['1KL 


o 
o 

Y 


1        0 
0        cos  Jj 
0       -sin  $ 


0 

' 

• 

+ 

n 

0 

0 

-sin  6 

cos  G  sin  j> 

cos  6 

cos  ^ 

=  T, 


II  -  3  Conservation  of  Momentum  Equations 

The  derivation  of  the  dynamical  equations  for  a  vehicle  with  a 
varying  mass  and  center  of  mass  have  been  treated  in  two  different  manners 
by: 

(1)  Albert  Strumpf  of  Davidson  lab,  (see  ref  2)  using  vector 
calculus,  and  including  variations  in  mass,  moments  of  in- 
ertia and  CG  position, 

(2)  Pierre  Dog an  of  MIT  Instrumentation  Lab,  (see  ref  5).  using 
a  Langrangian  formalism,  includii  g  1  *  j  variation  of  the 
inertia!  tensor,  and  making  assumptions  as  to  the  form  of  the 
movable  weights. 

The  treatment  by  Strumpf  assumes  a  body  that  is  a  rocket,  so  that  mass  is 
discharged  from  the  body.  The  treatment  by  Dogan,  for  the  Deep  Submergence 
Rescue  Vehicle  (DSRV),  assumes  that  th  \     ss  of  the  body  is  constant  but 
allows  the  position  of  t)   -   '   i      ity  to  change. 

The  development  b,  ;      voids  t:   .'    \   I    ":  c   '   tions 
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involved  in  Strumpf's  development,  while  sacrifa  Hty  by  not 
incli  '  ,  a  change  in  mass.  This  feature  could  be  included  if  it  were 
desired  and  so  the  approach  by  Dog;  ,  \  ill  be  used. 

Nine  generalized  velocities  and  coordinates  are  sufficient  to 
describe  the  motions  of  the  vehicle.  These  are:  six  velocities  (u,  v,  w, 
p,  q,  r)  to  describe  the  vehicle  and  three  coordinates  (x>,  y  ,  zQ)   to 
describe  the  motion  of  the  center  of  gravity.  The  following  six  Lagrange 
equations  will  give  the  needed  vehicle  momentum  equations.  (See  Lamb, 
"Hydrodynamics*1  6th  Ed,  page  168) 

A^!y„  r4-;-q4  =  x  (20 -1) 

dt  -au     bv   4  "9W 

d  Hx  „  p"?Iv    ?Jv   y  (2,3  .  2) 

d  L  'civ     ^  3\:   '  J   'iU     j  ^ 

d  "*TV  7>T  -JT  , 

SEW  -  q  TiT+  PTvB  z 

d   3T„  ^T  "3T  ^T  "3T  v 

on-  3p  nv  ?w  "aq     .        ^r 

d  ^Ty  ^Tv  >Tv  ?vv  ?v.  2 

d    2>T„  2T,r  sT  2TV  2T 

s  -s?  -  v  in?  + « "^  ■•  q  "5F  +  p  9?  - K  (2-5  - 6) 

T-  is  the  vehicle  total  kinetic  energy.  X,  Y,  Z,  K,  M,  N  are' 
the  generalized  forces  and  moments  of  which  some  can  be  further  described 
by  a  potential  function  and  others  represent  friction  and  drag( 

The  three  Lagrange  equations  necessary  to  describe  the  momenti 
balance  equation  for  the  jettisonable  ballast  subsystem  are  assumed  to 
reduce  to  quasi- steady  equations  defining  the  CG  location  from  the  integrals 
of  the  various  ballast  release  rates. 

This  proceedure  is  as  much     ossity  as  it  is  a  simplific 
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since,  in  the  case  of  dropping  ballast,  tho  ballast  may  take  almost  any 
form  from  liquid  mercury  and  iron  shot  to  large  blocks  of  metal  or  pieces 
of  equipment •  The  dropping  of  a  liquid  or  a  granular  solid  can  be  reason- 
ably modeled  as  a  function  of  time  but  the  dropping  of  chunks  of  metal  or 
pieces  of  equipment  would  create  singularities  in  a  functional  relation- 
ship.  The  obvious  answer  would  be  to  use  a  combination  of  a  smooth  func- 
tion and  steps  to  obrain  a  reasonably  accurate  model  of  the  deballasting  of 
a  vehicle* 

There  is  one  additional  factor  which  also  affects  the  decision  to 
reduce  the  function  to  a  quasi-steady  process*  The  total  ballast  dropped 
is  no  more  than  three  percent  of  the  total  vehicle  weight  and  tho  rate  at 
which  it  is  removed  is  of  order  »3  percent  of  the  total  vehicle  weight  per 
second.  This  would  then  say  that  any  Taylor  series  expansion  of  this 
function,  which  retained  terras  commesurate  with  the  second  order  expansion 
to  be  used  in  obtaining  the  hydrodynamic  forces,  would  contain  at  most  the 
linear  terms. 

The  effect  of  such  an  approximation  is  entirely  dependent  upon 
the  length  of  the  time  interval  over  which  the  process  is  assumed  to  be 
steady  and  shall  be  discussed  in  conjunction  with  the  computer  program. 
Suffice  it  to  say  here,  that  until  accurate  model  tests  can  be  conducted, 
the  effect  of  this  assumption  is  truly  unknown,  but  appears  to  be  of  incon- 
sequential magnitude , 

The  equations  describing  the  variable  position  of  the  CG  in  the 
body  axis  system  are: 
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x   W 
GB 


K    t 


*b 


lsr. 


(2.3  -  ?) 


N    t 

yGBw-2y,  S  ^  « 


y^_  = 


i 


i  r 


W  -  >.' 


(2.3  -  8) 


%BW'-^«i|  Widt 


z,-,  = 


w 


2w 


(2.3  -  9) 


whore  x__,  y,.,-,t  z   are  the  components  of  the  CG  of  the  vehicle  with  all  of 
GB   GB   GB 

the  ballast,  W  is  the  vehicle  weight  including  all  jettisonable  ballast, 

x. ,  y. ,  z     are  the  CG's  of  the  N  jettisonable  ballast  weights  W  , 
i   a   i  i 

The  vehicle  kinetic  energy  is  the  sum: 


rn  -*o      — » 

T  =7t  +  fceol  M 
V        c.      G  G 


(2.3  -  10; 


where  n,  V  ,  I  and  co  are  the  vehicle  mass,  the  CG  velocity,  the  inertial 

G    b 

tensor  about  the  CG  and  the  angular  velocity  vector.  Defining  exes  x' ,  y' , 
£*  through  the  CG  parallel  to  the  vehicle  body  axes,  one  has 


JG  = 


""Vx* 

I 

X»y* 

I 

Vx' 

Vy' 

I 

I  . 

a'x' 

aVy' 

I 

'v.' 


I  •• 


»  »  • 


(2.3  -  11) 


■  v  re  -£  i    ft  I   ,    ,    ...,e  have  the  usu  laningsi 

xx        y  j' 


1'! 


(y*  "  +   z9    )   c^ 

v 

Coi..^  itirig  the  components  of  V     in  body  axes  one  has! 

v 

VQ  =  (u  +  qzG  -  ryQ  •:-  %,  v  +  rxG  -  pzG  +  yQ,  w  +  pyQ  -  qx£  +  zQ) 

(2.3  -  12) 

Substituting  (2,3  -  11)  end  (2.3  -  12)  into  (2,3  -  10)   one  gets 
the  ld.netic  energy  T     and  its  partial  derivatives: 

_J£  =  m  (v.  +  qzQ  -  ryQ  "+  xQ)  (2.3  -  13) 

^-"(v+nfe-P^  +  V  (2.3-14) 

^  =  m  («  +  PyG  -  qzG  +  y  (2.3  -  15) 

•—---  =  1  p  •:-  I  q  +  J    ,     -r 

"<?i  x'x'  x*y  x'z' 

-  mzQ  (v  •!-  rxQ  -  pz    +  »yQ) 

+  myG   (v  +  pyQ  -  qx    +  ZQ)  (2,3  -  16) 


3TV 

=  I    ,    ,q  +  I    ,    ir  +  I    ,    ,p 

"aq        yy  y  »  y'x* 

"  r^  (w  +  wrG  -  <*XG  +  *g> 

+  mzQ  (u  +  qzG  -  ryQ  +  xj  (2,3  -  1?) 

+  inxG  (v  ■:   r  ■    ~  pz    +  yQ)  (2,3  ~  18) 
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Tho  inertial  tensor  I-  about  the  CG  can  bo  algel     lly  related  to  tl 

G 

inertial  tone or  I  about  th<  incd  in  x,  y,  z)  and  the  disk 

o 

X .  i  y  i  z  •  For  example : 
G   G   G 

J   =  I  .  ,  •!-  m   (s2  +  v2) 
-xx    x'x'     V"G   V 

1=1     -  mx„y„ 
xy    x'y*     G  G 

Substituting  these  for  the  inertial  tensor,  1„ ,  equations  (2,3  -  16) 
through  (2.3  -  18)  becomo: 

-—  =  I  p  •!-  I  q  +  I  r 

^)P    xx^    xy1    xz 

-  mzQ  (v  +  y  ) 

+  mv,  (w  +  z„)  (2,3-19) 

b        u 
— 2  =  1  q  +  I  r  +  I  p 

^q.   yy   yz   jrc 

-  mx_  (w  -!-  sL) 

+  mz_  U  •:  x  )  (2.3  -  20) 

— X  =  I  r  +  I  p  *  I  q 
e>r    fcz     zx     ZJ 

-  myG  (u  +  xQ) 

+  mx,  (v  +  y  )  (2.3  -  21) 

Substituting  the so  equations  and  equations  (2.3  -  13)  through  (2,3  -  15 ) 
(irithout  primed  subscripts)  into  the  Lagrangian  equations  of  motion 
(2o3  -  1)  through  (2.3  -  6),  one  gets  tho  dynamical  equations: 
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X  s  m  [u  ~  rv  +  qw  -  x     (q2  +  r2)  +  yQ   Cpq  ■     ' )  *    zQ  (pr  +  4) 

+  2qSr   -  2rv,   +  xJ  (2-3-  22) 

^  G  G         Gj 

I  =  m  [v  -  pw  +  ru  -  yQ  (r2  +  p2)  +  zQ  (qr  -  p)  +  xQ   (qp  +  r) 

+  2rxQ  -  2P5G  *  yc]  (2.3-23) 

Z  =  1.1  Pw  -  qu  +  pv  -  z     (p2  +  q2)  +  xQ   (rp  -  q)  +  yQ  (rq  +  p) 


+  ?PyG  -  2qV+  V) 


(2.3  -  #0 


r  (q    ■  Pr)  +  I_  (r  'i-  pq)  '  (q2  -  A  +  (*sz  -  *yy>  qr 


xx1  xy 


xz 


yz 


-  mz,,   (v  -  pw  +  ru  +  ix,  +  yr) 

b  Ct  b 


+  my     (w  -  qu  +  pv  -  qx     +  \) 


xx  xy 


}  ■'.  ■ 


•  •  9 

•S-  I     p  +  I     q  +  I     r  +  rax,..   (z  r  +  y  a) 

G        G  G 

M  =  I     q"  +  I        (r  -  qp)  +  I        (p  +  qr)  +  I    r  (r 


yy         yz  yx 

-  roxu   (w  -  qu  +  pv  +  py    +  z   ) 

■'    i    .     ( '    •    j  \    ■!    q.;   -•   :•.'..   ••  '    . ) 


zx 


(2.3  -  25) 


p^)  +  (I       -  I     )  rp 

r  xx         zz 


G 


G       -'G' 


-I-  1     q  +  I     r  +  I     p  +  i.iyp   (x-p  +  z  r) 


yy        y^        yx- 

N  =  I_r  ■:-  Izx  (p.-  rq)  +  IRV  (q  'I-  rP)  +  1^  (p 


2 


zs 


sy 


(2.3  -  26) 


12)  •'    (Iw  -  W   ^ 

yy       x*- 


my^   (u  -  rv  +  qw  +  qz.     -I-  x.. ) 

b  b  b 


+  rnx^  (v  -  pw  +  ru  -  pz     +  y  ) 

+  I     r  +  1     p  +  I     a  +  mz     (t  a  +  x  n) 

zz  zx^         sy1  G     JQq       ^gP 


(2.3  ~  2?) 


1? 


wh  re  I  ,1   . ....  are  about  axes  through  the  vehie"1   origin*  T] 
xx   Sy  t 

equations  also  contain  derivatives  of  tho  form  I  ,1   **•••  XL.  x  ••••••. 

a  xx   xy       C :   G 

The  tensor  I  can  be  thought  c     being  made  up  of  two  parts;  .• 
o 

constant  part  representing  the  vehicle  with  all  jettisonable  ballast  re~ 
moved,  I  ,  and  the  jettisonable  ballast,  I  ,  such  that: 

I  =  1+1,.  (2.3  -  28) 

o    1    ^ 

The  contribution  of  the  ballast  can  be  modeled  by  assuming  that  each  weij 
that  is  part  of  the  jettisonable  ballast  is  lumped  at  its  center  of  gravity 

Cx  . ,  y  .  z  , ) . 

gi   gx   gn 

Then: 

I  "    =  3  ri--  ^ •  '''•  sZ.)  (2.3  -  29) 

xx,2   5  ^  i   gi    gi 

N 

1     =-  >   m.x  .y  (2.3  -  30) 

xy,2   iTl  l  S:'-  gi 

N 

I   ^  «  >   m.x  .z   .  ...   '  (2e3  -  31) 

xz,2   i  =  i  1  e1  E1 


where  N  is  the  number  of  weights. 

1  is  the  time  varying  part  of  I  ,  however,  to  be  consistent  wit] 
t]   quasi-steady  approximation  made  in  developing  equations  (2.3  -  7 
through  2*3  -  9)»  a  quasi-steady  change  in  the  inertial  tensor  must  also 
be  assumed.  This  then  says  that  the  terms  involving  ti     ;  derivatives  of 
the  inertial  tensor  and  the  center  of  gravity  can  be  dropped,  since  they 
are  zero  during  the  time  interval  over  which  the  process  is  assumed  to  be 

The  equations  can  be  furth     3    3  when  it  i.r;  : 
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duo  to  the  assumption  of  an  xz  piano  of  symmetry  and  a  fu:  '   r  t  tion 

that  the  balla*  b  will  be  dropped  symmetrically,  the  products  of  inert:! 

I       :  T     .  B  I     ^  =  1       :  1       al       =  I       =  I       =  0 

xy,l    yx,l    yz,l    zy,l    xy{2    yx,2    yz,2    zy,2 

Thus,  the  total  moment  of  inertia  becomes  quasi- steady  and  can  be  used  in 
the  equations  of  motion  -without  0,  i  or  2  subscripts.  The  final  form  of 
the  dynamical  equations  then  becomes s 

X  s=  raf  u  -  rv  +  qw  -  x  (q  +  r  )  +  yG  Cpq  -  r)  +  Zn   Cpr.+  q)i 

(20  -  32) 

Y  =  m  [v  -  pw  +  ru  -  yQ  Cr2  +  p2)  +  ^   (qr  -  p)  +  x^  (qp  +  r)] 

(2.3  -  33) 

Z  =  n  [w  -  qu  ■!•  pv  -  zQ  (p2  +  q2)  +  x  (rp  -  q)  +  yQ   (rq  +  p)j 

(2,3  -  3;i) 

K  =  I  p  +  I   (r  +  pq)  +  (I   -  I  )  qr 

xx    xz  2-21   yy 

+  m  ["-  ^G  (v  -  pw  +  ru)  +  yQ  (w  +  pv  -  quj]       (2,3  ~  35) 

M  =  I  4-1-1   (r2  -  p2)  -1-  (I   -  I  )  rp 
yy    ex  xx   zz 

+  m  [--  xG  (w  -  qu  +  pv)  +  zQ  (u  +  qw  -  rv)]       (2,3  -  36) 

N  =  J  *  +  I   (p  .  rq)  +  (I   „.  I  )  ^ 

zz    zx      *     yy    xx   ^ 

+  m  T~y  (u  -  rv+  qw)  H  :  ^    (\  +  ru  -  pw)J        (2.3-3?) 
where  m,  x  ,  y  ,  z  .  I  ,1  ,  I  ,  I   are  all  quasi-steady  functions  of 

u    U    U    XX    yy    zz    X'/, 

time , 

The  forcing  terms  for  the  hydrodynamic  equation:  ■     3  up  of 


19 


gravity  forces,  hydrodynamic  forces  i  1  propulsion  induced  forces.  Since 
the  problem  has  been  defined  as  a  free  ascent  pi      £  trill  be  no 

propulsion  forces  included.  Because  of  this  restriction,  the  expe. 
hydrodynamic  coefficients  should  be  obt:  '    without  propellers  running. 
The  total  forcing  i  are  then: 

X  "  XG  "'   XH 
Is=YG  +  YH 

z  .  zG  +  zH 

N  =  NQ  +  NH  (2.3  -  38) 

^  "  ^  Gravity  Forces 

The  hydrostatic  pressure  field  induced  by  gravity  creates -a 
buoyancy  force  B  through  the  CB,  This  force  varies  with  ambient  water 
density  and  vehicle  volume.  The  instantaneous  weight  of  the  body,  W,  act- 
ing through  the  CG,  is  made  up  of  W  ,  the  weight  of  the  body  trithout  jetti- 

sonable  ballast,  less  W  ,  the  weight  of  the  ballast  components  removed. 

i  N 

During  the  initial  phase  of  the  ascent  the  items  in  ^>   '  VJ. 

i  =  i 
rro  increased  until  all  of  the  ballast  components  have  been  removed,  at 

this  point  the  buoyant  force  (-B  +  W)  becomes  a  maximum,  if  we  neglect 
changes  in  vehicle  volume  and  density.  This  maximum  force  is  sustained  for 
the  rest  of  the  ascent. 

The  instantaneous  buoyant  forco  can  thus  be  represented  by: 
N 
i  ~   1  ~ 
^v    1  :  the  nu  ]        of  ballast  c<     nts  removed.  This  force  acts 
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upwards  along  tho  local  vertical.  Da©  to      sice   of  origin  oi       dy 

s,  the  gravity  induced  torqu  ii  r  :  B  where  r  is  the  CG  position 
vector  in  body  axes  and  B  is  the  vehicle  buoyancy  vector  acting  up  along 
the  local  vertical,,  Resolving  along  body  axes  tho  gravity  forces  beco:;:o: 

T     =  -  (W  -  B)  sir.  0  (2  .4  -  2) 

G 

Y  ^  (u  -  B)  cos  C  sin  $  (2.4  -  3) 

o 

ZQ  -  (w  -  B)  cos  0  cos  j>  (2.4  -  4) 

KQ  =  (y  W  -  ypB)  cos  G  cos  <f>   -  UQ*W  -  z^B)   cos  C  sin  j>   (2.4  -  5) 
M  =  ~  (xlW  -  xJB)  cos  0  cos  o7  -  (zrW  -  z  B)  sin  0     (2.4  -  6) 

G        Tj      B  'It      h 

NG  =  (xJJ  -XgB)  cos  G  sin  0  +  (yGW  -  yBB)  sin  G      (2.4  -  ?) 

II  -  5  Hydrpdynamic  Forces. 

The  hydrodynamic  forces  and  moments  that  act  on  a  body  moving 
through  a  real  fluid  are  the  result  of: 

(1)  hydrodynamic  inertial  affects  (linear  added  mass  terms), 

(2)  skin  friction,  separation  and  cross-flow  drag  effects, 

(3)  circulation  effects. 

These  effects  are  all  functions  of  the  velocities  and  accelerations  of  tho 
body.  Therefore,  the  hydrodynamic  forces  and  moments  can  be  expressed 
f ui  ic tionally  as : 


-=  f  (u,  v,  irt  p,  q,  r,  u,  v,  w,  p",  q,  rj  (2.5  ~  l) 


This  function  may  bo  reduced  to  a  workable  form  b;>.  <  ]  mding  the 
function  in  a  Taylor  series.  Expanding  the  functic  in  this  form  requires 
that  the  f unctii   and  i1  '  nuous  in  tl 


values  of  the  variables  under  consideration.  A  typical  i      order  t    . 
sion  of  one  of  the  force  equations  would  bo  of  the  form: 

v  _  y  ..  „  ^  j.   *^X  .  ■   "2 

A    =    A         t     U    — — -     '      V   — *  •-     '       tceeestseecococcc  •     J-     - 

O  -J,U  5>V  3r 

,      •   *X     ,      .    >X     .  ...     •    ?-' 

2>u  *v  *r 

.    x    /   2  A    ,      2   >?:     ,  .    .2  A 

T    V     \U       — — iy     >     V        '    '    _     *'       iitcettiic  *'      I'       ~~~~ *y 

^u  "b  V"  t>  r 

2  ? 

"bX     J  ^X       .  ,       ,  ^j2X 

+  uv +  uw —  +   .  e  c  o  e  c . .  •  t  o  <■  ui*  -~~r-~ 


+  pq_  +  pr_  +  qr_ 

+  uu  —  — b-  +  uv  —rrx  +   * « <  •  c  <  »  i  *  t »  +  ur  — -—v 

}U}U  SU^^  3UiV 


•  *x    ,  ,    .  ^ x 

'        -'    l  '  C         •         0    C    i    C    C    C    t   C   C    C    C    C    t    «    C   C    0    4   i    C    £  *JLX  • 

r  u  ^r^r 

• .   "*  X  . ,    5  X  . .    -3  / 

+    UU   ~2"    +    W   ""72    +     •  •  •  c  .  e  .  t  e  .  .  *  +    IT    •    v 


t      •  •     cl     X         (  •  t     i>    A  x  •» 

*•"     PQ  •        •"     *1        •*Ct««(CCi«C*O*0«0**«      "*"     Qr  ■        *  li-l^      *"      *i  / 

ap^q  ^r 

This  equation  contains  ninety-one  constant,  linear  and  second  order  terras 
arising  from  the  Taylor  series  expansion.  The  number  of  terms,  however, 

be  reduced  to  thirty-three  on  the  basis  of  the  problem  assumptions  and 
restrictions  delineated  in  Chapter  I, 

The  constant  term  is  d  ]   3  to  confoi   Lth  ment  that 
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initially  the  only  disturbs-,   :'■  c   to  a  buoyant  force  which  is  included 
as  a  gravity  term. 

The  terras  involving  products  of  accelerations  viith  velocities  or 
accelerations  are  dropped  because  the  second  order  effects  vrere  restrict 
to  velocities  only.  In  addition  the  results  of  potential  theory  indicate 
that  these  derivatives  are  aero,   (See  ref  2) 

After  el:'  i   ting  all  but  the  linear  terms  and  the  second  order 
velocity  terms,  the  Taylor  series  expansions  of  the  hydrodynamic  forces  and 
i   ants  become  of  the  form: 


X  =  X  u  +  X  v  +  X  w  +  X  p  -:-  X  q  +  X  r 
H    u     v    w     p^    q^    r 

+  -X,u  +  X.v  +  X.w  +  X.p  •!•  X.q  +  X.r 
u     v    w.    p     q     r 

+  X  2u2  •;-  X  2v2  +  X  ^   +  X  p2  +  X  a2  +  X  2r2 
u      v      vt      p      q~     r 

+  X  uv  +  X  "uw  +  X  up  +  X  uq  +  X  ur 
uv     uw     up     uq      ur 


+  X  v\r  +  X  vp  +  X  vq  -s-  a  'v 

vw       vr»       vn 


-5-  X  '• 

vw     vp  *  vq      vr 


+  X  wp  +  X  v?q  +  X  wr 
wp  c        vrq  ^    wr 

•{-  X  pq  +  X  pr  +  X  qr  (2.5  -  3) 

Pq      pr      qr 

The  terms  in  this  expression  are  seen  to  fall  into  one  of  three 
catagories,  namely: 

(1)  added  mass  or  inertial, 

(2)  second  order  non- inertial, 

(3)  linear e 

In  order  to  f  1  .  :  ]         the  num.'    of  1       '  'ned  in    h  < 
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'  >,  it  is  nee  sary  to  look  at  the   '       ;:;o  terms  and  the  effc 
that  a  p? .   of  symmetry  has  on  them. 

II  -  5«i  Added  Mass-  Terms 

A  body  moving  through  a  real  fluid  induces  a  motion  in  the  other- 
wise stationary  fluid  because  the  fluid  must  move  aside  and  then  close  in 
behind  the  body.  As  a  result  of  this  motion  the  fluid  possesses  kinetic 
energy  that  it  would  not  possess  if  the  body  were  not  in  motion.  The  added 
mass  terms  in  the  equations  take  into  account  the  energy  given  to  the  fluid 
by  the  body. 

If  the  body  motion  is  steady,  the  related  fluid  motion  is  also 
steady  which  requires  that  the  kinetic  energy  be  constant.  If  the  kinetic 
energy  is  constant,  no  work-  is  being  done  on  the  fluid  and  therefore  the 
added  mass  terms  may  be  emitted. 

If,  however,  the  body  is  in  accelerated  motion,  there  will  be 
work  done  by  the  body  on  the  fluid  and  it  will  be  necessary  to  retain  at 
least  some  of  the  added  mass  terms. 

Work  is  accomplished  by  moving  a  force  through  a  distance.  In 
the  case  of  a  submerged  body,  the  distance  is  the  distance  the  body  travels 
and  the  force  is  the  integral  over  the  surface  of  the  body  of  the  pressures 
exerted  by  the  body  on  the  fluid.  This  force,  in  general,  represents  a 
system  of  forces  and  moments  acting  on  the  body  which  can  be  obtained  from 
equations  (2.3  -  i  through  6),  when  the  kinetic  energy  is  varied. 

The  kinetic  energy  of  the  fluid  can  bo  represented  as  a  function 
of  the  six  velocity  components  (u,  v,  w,  p,  q,  r) .  A  quadratic  form  of  this 
function  as  given  by  Lamb  (see  page  1?2  ref  i)  is: 


: 


2T  =  Au2  •:-  Bv2  +  Cw2  +  2A'  <n  ■:  2B'wu  +  2C*uv 
+  Pp2  -1-  Qc.2  +  IV-2  +  2P'qr  +  2Q'rp  +  2R:; 
+  2Lup  +  2IIvq  •!•  2'  ' 

•:-  2F  (vr  +  wq)  +  2G  (wp  +  ur)  +  2H  (uq  ■:■ 
-i-  2F'    (vr  -  wq)  +  2G*    (wp  ~  ur)  +  2H'    (uq  -  vp)      (2.5  -  *0 

where  the  twenty-one  coefficients  A,  Bt  C  etc,  are  certain  constants  deter- 
mined by  the  form  and  position  of  the  surface  relative  to  the  co-ordinate 
axes. 

Letting 


F 

+ 

F' 

ri 

F 
1 

H 

+ 

II  • 

= 

Hi 

G 

+ 

G' 

= 

Gi 

F  ~  F' 

=  F2 

H  -  II1 

mBz 

G  -  G* 

z,  G 
2 

the  last  six  terms  in  (2,5  -  ^)  becomes 
+  2F  vr  +  G  sp  •!•  21-^  uq 
+  2F2wq  +  G2ur  +  2H  vp 

Nine  of  the  coefficients  in  the  above  expression  may  be  sot  to 
zero  if  we  take  advantage  of  the  assumed  xz  plane  of  symmetry.  Symmetry 
arguments  say  that  T  should  not  be  changed  if  any  of  the  terms  vw,  uv,  qr, 
pqf  up,  vq,  vr,  wp,  ur  is  replaced  respectively  by  (-v)  w,  (-v)  u,  (-r)  q, 
(~p)  qi  (-p)  u,  (~v)  q,  (-r)  w,  (-p)  w,  (-r)  u.  Those  terms  correspond  re- 
spectively with  the  coefficients  A',  C,  P',  R*t  L,  m'  n'  G,  G«. 

When  T  is  substituted  for  T  in  equations  (2,3  -  1  through  6), 
these  nine  coefficient  d  to.  their  corresponding  '     in  equa- 

ls (2»5  •■  3)(  at  which  ti        t 
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hydr<  i;        tic  expressions  0 

Six  partial  derivati .  '    be  obtained  J  '  Ion  (2.5  -  ^) 

in  order  to  expand  equations  (2.3  -  1   through  6),    these  are:. 

t±L  =  Au  +  B'w  +  Cfv  +  Lp  +  11  q  +  G  r  (2.5  -  5) 

~-£  s  Bv  +  A'tr  +  C«u'+  Kq  +  F  r  +  H  p  (2.5  -  6) 

<J  V  X  & 

—  =  Cw  +  A'v  +  B'u  +  Nr  +  G.p  +  F0q  (2.5  -  7) 

■^  =  Pp  +  Q'r  +  R'q  +  La  +  G^j  +  H  v  (2.5  -  8) 

— -  =  Qq  -i-  P'r  +  li'p  +  Hv  +  II  u  +  F„w  (2.5  -  9) 

"^q  12 

~d  T 

— £  =  Kr  +  P*q  +  Q'p  +  Nw  +  F«  v  +  G„u  (2.5  -  10) 

^  r  12 

Those  derivatives  are  then  substituted  in  equations   (2.3  ~  1  through  6)  to 
obtain  : 

X  =  Au  +  B'w  +  C'v  +  Lp  +  H,q  +  G  r 

iH         2 

+  Cqv  +  A'vq  -I-  Bcuq  +  Nrq  -I-  G  pq  +  F  qq 

-  Bvr  -  A'wr  -  C'ur  -  Ilqr  ~  F„rr  -  H  pr  (2.5  -  11 ) 

1  2 

I  =  Bv  +  A'w  +  C'u  +  llq  +  F  r  +  H  p 

l  ^ 

+  Aur  +  B'wr  J<-  C'vr  +  Lpr  +  H  qr  +  G_rr 

-  Cwp  -  A'vp  -~  B'up  -  Krp  -  G  pp  ~  F  qp  (2.5  -  12) 
Z  rn  Cv  +  A'v  +  B'u  +  N*  +  Gp  -I-  F  4 

+  Bvp  +  A'wp  +  C'up  -I-  Mqp  +  F  'rp  ■+  H  pp 

-  Auq  -  BVq  -  C'vq  -  Lpq  -  H  qq  -  G^,  q  (2.5  -  13) 
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K  a  Pp  -:-  Q'r  +  R'q  +  Lu  +  G,w  +  1  ./ 

•!■  Rrq  H    P'qq  +  Q'pq  H    Nwq  +  F  vq  +  G  uq 

-  Qqr  -  P'rr  -  R'pr  -  Mvr  -  H  ur  -  Fgi 

•{-  Cirv  +  A'vv  +  B'uv  -l-  Nrv  +  G  pv  +  F?qv 

-  B       -  A'ww  -  Ceix;:-  -  Mow  -  F„  to  ••  II  pw  (2.5  -  Vl) 

\.  -  1  2 

M  =  Qq  +  P':    •:   R'p  +  Mv  +  H  u  +  F^t 

+  Ppr  -J-  Q'rr  +  R'qr  +  Lur  +  G  wr  +  H  vr 

-  Rrp  -  P'qp  -  Q'pp  -  Nvip  -  F  vp  -  G  up 
+  Auw  +  B'ww  +  C'vw  +  Lpw  +  H/qw  +  G  n: 

-  Cwu  -  A' vi.  -  B'uu  -  Kru  -  G  pu  -  F  qu  (2,5  -  15 ) 

N  =  Rr  -J-  P*q  +  Nw  +  F.v  +  G-  u 
1  1  2 

+  Qqp  +  P»rp  +  R'pp  •:   Myp  +  Hup  +  F  wp 

-  Ppq  -  Q'rq  -  R'qq  -  Luq  -  G  wq  -  H.vq 
+  Bvu  +  A'wu  +  C'uu  +  Mqu  -i-  F  ru  +  H  pu 

J.  Cs 

-  Auv  -  B»wv  -  C'w  -  LpTj  •  H^qv  -  G  rv  (?,«5  -  16 ) 

Setting  the  various  coefficients  in  equations  (2,5  ~  H  through 
16)  equal  to  their  counterparts  in  equations  (205  -  3)  it  is  found  that? 


2? 


A'  a  X 


vq 


X    =  Y.  a  -  Y    a  Z.  a  Z    =  K    =  -  K 

wr    w      vp    v    vp    w      wi      vu    "wu. 


C  a  X.  =  -  X    a  I.  =  3     =  Z    =  -  Z     :   -   I  =  N     =  -  N 

v      ur    u    vt    up      vq      uw    vw    uu      w 


S^pq^ 


f   a  Z,  =  K.  =  K   =  M 

PP     1     w     pv     X  . 


=  ~  M 


^  ~  N 
1  '      wq 


G2  "  Xr  =  Yrr  = 


Z   =  K 
rq    uq 


vv    rw    u 


5  v 


L  a  X.  a  K.  a  3    a  H   a  -  K   a  -  1! 
p    u    ur    pw     uq     pv 


M  a  Y.  =  -  K    a  M,  a  N  .   a  K 

q      -qw    v    vp    qu 


N  a  Z. 


K    a  -  M    a 
wq      wp 


K     a  NB 

ru    w 


P»  =  K   a  -  K   =  M.  =  M   =  N,  r  ]; 

qq      r,     r    qr    q    rp 


R»  a  K.  a  ~  K    a  M.  B  -  M    =  K    a  _  K 

q      pr    p      qp    pp      qq 
/    / 

K  -  K   a  X    a  K 

qr    vr 


L  -  N  =  Y   =  M   b  M 

pr    ru    pw 


M  -  L  =  Z 


pq 


G  -:-  go  =     a  -  K 

i    2    j,'      pu 


li  -  L  a  N     =  N 

vp    uq 


(2.5  -  1?) 


H,  +  H   a  13     -.N 

1    2    up      Vq 


B  -  A  =  N 


V; 


(2.5  -  18) 
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A    =  X.    =   Y  :      -    Z 

u         ur  uq 


B  =  -  X       =  Y.  =  Z 

vr         v         vp 


C   =  X 


qv 


Y       =  Z. 

vp        w 


B'  =  X.  =  X       =  Y 

vr         uq         vr 


Y       =  Z. 
up        u 


Z  r=    K  b] 

vq         uv         w 


M       b  -  N 
uu  ^  ■ 


r  b  i:,  -  »  i; 

P  pq 


Q  =  M.  =  N 

q      qp 


R  :=  N. 


Q»      B     K.      B     K  B    M  K  —     M  B     J 

r         )  -j         .1 ,  pp         p 


ro 


?*■  b  -  X       b  Y.  b  Z       =  -  M       =  N.  b  K 


rr 


rp 


vp        v        r-u 


F0  b  X       =  -  Y       b  Z.  b  M.  =  »  M       b  N 

2         qq  qp         q         VT  qu         vp 


II     =  X.   b  Y       b  -  Z 

1  q         qr  qq 


K       b  M.  s  H       =  N       =  -  N 


v. 


V. 


qv         up 


qv 


H     b  -  X       b  Yc  =  Z       b  K.  b  -  K       ■■    y       =  -  M       b  N 
2  pr         p         j,--,         v  p;-7         vr  vq         Pu 


C  -  B  b  K 


w 


F.   +  F     b  -  K       b  K 
1  2  vrr  vq 


R  -  Q  =  K 


qr 


P  -  R  b  M 


pr 


A  -  C   :    M 


vu 


Q    -   P    : 


1 


(2.5  -  18)  cont. 
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.•e  equations  (2.5  -  17)  represent  the  terms  whicl  to 

symmetry. 

Eliminj  •  '■  .,  those  terms  which,  due  to  the  xz  plane  of  symmetry, 
ai-o  zero,  and  using  tl  standard  nomenclature,  as  given  in  reference  (6), 
in  place  of  Lamb's  notation,  the  added  mass  expressions  becc: 


X.T.  =  X.u  •!■  X.w  -i\X,q  •!-  X  qw  +  X  uq 
Hi    u     w     q     qw      uq 

+  X  qq  +  X  vr  +  X  0r2  +  X  pi-  (2.5  -  19) 

qq     vr     j/     pr 

Y   =  I.v  •!  Y.r  +  Y.p  +  Y  ur  +  Y  wp 
Hi    v    r    p    ur     wr 

+  Y  qr  +  Y  wp  +  Y  up  +  I  qp  (2.5  -  20) 

qr     wp  ^    up     qpM^ 

Z  T.  =  Z.u  +  Z.u  +  Z.o  ■+  Z  vp  -!-  Z  r-p 
*«■    w    u     q     vp     rp 

+  Z  pp  +  Z  uq  +  Z  wq  +  Z  qq  (2.5-21) 

PP      uq      wq      qq 

K.  =  K.p  +  K.r  +  K.v  +  K  rq  -1-  K  pq 
Hi    P1     r     V     rq  x  pq^1 

+  K  vq  +  K  uv  +  K  ur  +  K  wr 
vq      uv     ur     \tt 

+  K  wv  +  K  pw  (2.5  -  22) 

wv     pn/ 

M«.   =  M.q  +  iicu  +  M.w  +  M     rr  +  M     vr 
ni        cl  u  v  rr  vr 

+  M     rp  +  M     pp  +  21     Vp  +  M     uw 
rp  pp1^        vp  uw 

+  M  ww  +  H  qw  -I-  M  uu  +  M  qu  (2.5-23) 

ms  qxr      uu      qu L  *u  •" 

N„.  =  K.r  +  K.p  +  K.-C  •:  K  qp  +  K  up 
Hi    r     p     v     QP      up  1 

+  K  .  wp  +  K  rq  +  N  vq  ■:  I   vu 
i  p  x    rq  J    vq  ^    vu 

+  Nrui;'  ':'  NwvMV  <2<: 
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The  equalities  given  in  <      ns(2#5-l?  and  18)  are  based 
entirely  upon  potential  theory  and  do  not  necessarily  hold  in  the  presence 
of  circulation  and  viscous  effects.  If  circulation  and  viscous  effects  are 
neglected,  these  equalities  provide  reasonably  good  estimates  of  some  of 
the  second  order  coefficients  which  are  not  directly  amenable  to  measure- 
ment by  conventional  towing  tank  techniques, 

II  -  5*2  Second  Order JT^erras 

The  second  order  terms  which  are  no  I  ■   Ided  mass  terms  are: 

XH2  =  X  2u2  +  X  2V?'  +  X  2ir2  +  X  2p2 
u      v      v      p 

+  X  uv  *!-  X  uw  +  X  up 
uv      uw      up 

+  X  vw  +  X  vp  •:-  X  wp  (2.5  -  25; 

vw     vp     wp 

hz  -   Y  2»2  +  Y  7?   +  Y  /  +  Y  / 

u      v      w      q 

+  Y  uv  +  Y  uw  •!-  X  uq 

uv      uw      uq 

+  Y vw  +  Y  vq  +  Y  wq  (2.5  -  26) 

™  vq     wq 

ZH2  ■  Z  2»2  +  Z  2**  +  Z  /"  -1  Z  / 

u      v      w      r 

+  Z  uv  •!•  Z  uw  +  Z  ur 
uv     uw     ur 

+  Z  vw  +  Z  vr  +  Z  wr  (2.5  ~  2?) 

vw     vr     wr 

KH?  =  K  2u2  +  K?P2  +  K,,„UP  (2--5  "  28) 

^    u      p^     l1P 

J  i;2  "  M  2^  +  M  2^  +  \Q^  &'5   -  29) 

v       q       v^ 

%2=N2w2  +  M2r2H    \      '-      '  (2.5-30) 

v  1- 
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The  symmetry  plane  (xz  plane)  .force  and  ]  (X,  Z,  M), 

which  involve  products  of  the  syi   tr;  p]  ne  velocities  (u,  vr,  q)  Kith  the 
out  of  piano  velocities  (v,  p,  r)  rust  be  zero  for  the  reason  that  the  same 
force  o,  ;    nt  must  result  if  v,  p,  r  are  replaced  by  -v,  -p,  ~re  This 
does  not  hold  true  with  respect  to  products  of  xz   plane  velc.  '     or  to 
products  of  out  of  plane  velocities . 

Another  symmetry  t    ,    nt  that  can  be  made  is  that  .:  symmetry 
plane  velocity  or  products  of  symmetry  plane  velocities  should  not  cause 
out  of  the  plane  forces  or  moments.  For  example,  there  should  be  no  Y 
force  resulting  from  a  w  velocity  or  even  a  combination  of  w  and  u  ve- 
locities , 

After  removing  these  syj   I  ;  \      as  from  equations  (2»5  -  25 
gh  3°)  the  second  ord  r  hydrodynamic  terms  left  a] 

XH2  =  X  /  -■   X  /  +  X  /  +  X  2P2 

u      v      w      p 

+  \,^ +  V*     ,  <2-5  -  3D 

YH2  "  X/'  *   V  +  V"  +  V  '    (2"5  "  32) 

ZH2  "  Z  /  +  Z  2v2  '•  Z  /'  •'  Z  2rZ 
u      v      w      r 

•:•  Z  uw  +  Z  vr  (2„5  -  33) 

uw     vr 

KH2  ^  */"  +  VP  (2<5  "  3'} 

*H2  "  M  2v2  +  M  2^  (2*5  "  35) 

V         q 

NH2  =  N  2r2  +  VW  (2'5"  * 

r 
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II  -  5*3     Linear  Ten 

The  linoar  terms  in  the  hydi-odynamic  force  equations  are: 

X1TT   =  X  ii  +  X  v  +  X  w  +  X  p  +  X  q  +  X  r  (2.5  -  37) 

IIJ-*        u  v  ^  p  q  r 

yut  =  y  u  +  i  v  h  y  w  +  y  p  +  y  q  +  i  r  (2,5  -  38) 

riL        u  v  w  p  q  r 

Zur   =  Zu+Z  v  +  Z  w  +  Zp  +  Zq  +  Zr  (2,5-39) 

HI.        u  v  17  p^        qn         r 

K  T   s  K  u  +  K  v  -1-  K  w  +  Kp+Kq  +  Kr  (2.5  -  W 

HL         u  v  w  P  q  *" 

1L   =  M  u  +  M  v  +  M  w  +  M  p  ■:•  M  q  +  M  r  (2.5  -  **1 ) 

HL         u  v  w.  p  q  r 

Nu,   =  N  vi  +  K  v  +  N  v  +  N  p  -1-  K  q  +  N  r  (2.5  -  hZ) 

HL         u  v  w  p  q  p 

In  order  to  determine  which -of  the*   ter]   should  be  retained,  each  term 
must,  in  general  be  considered  on  its  own  merits.  There  is,  however,  one 
group  of  terms  all  of  which  may  be  eliminated  on  the  basis  of  the  assumed 
xz  piano  of  symmetry. 

The  coefficients  which,  due  to   symmetry,  must  be  zero  are: 

(1)  those  that  involve  derivatives  of  the  symmetry  plane  forces 
and  moment  with  respect  to  the  out  of  plane  velocities  (v,  p, 
r), 

(2)  the  out  of  plane  force  (Y)  and  moments  (K,  II)  which  would 
arise  from  an  in  the  symmetry  plane  velocity  (u,  w,  q). 

Those  that  fall  in  the  first  category  ar<    the  typo  which  re- 
quire the  force  or  moment  to  si  y  the  same  while  the  velocity  can  change 
sign.  Those  of  the  second  category  v,    ild  req  '     force  ]     licular  ' 
the  plane  of  s;    try  to  result  from  a  flow  in  the  plane  of  s;    try, 
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Er'  bing  those  terms  from  the  lin   .      te:        eq      '  Ions,  results 

in  the  following: 

XTvr   =  X  u  +  X  w  +  X  q  (2,5  -  h>3) 

HL         u  w  q 

Y      =  Y  v  +  Y  p  +  I  r  (2.5  -  W 

HL        v  p  r 


Zm   =  Z  u  +  Z  w  +  Z  q  (2.5  -  z:5) 

q 


JHL         u  w 


FHtt   =  K  v  +  K  p  +  K  r  (2.5  -  ^6  J 

HL        v  p^         r 

WU    =  M  u  +  M  w  +  M  q  (2.5  -  W 

tiL         u  w  q 

Nm   =  N  v  -:   N  p  +  N  r  (2.5  -  4-8) 

HJb        v  p  r 

II  -  5«^'  Discussion  of  Terns 

The  hydrodynamic  force  and  moment  equations  that  result  after 
eliminating  the  symmetry  terms  will  now  be  looked  at  in  order  to  further 
reduce  tho  number  of  terms  in  the  equations. 

Of  first  importance  is  the  decision  to  retain  all  terms,  or  a 
form  from  which  they  can  be  derived,  that  are  included  in  "The  Standard 
Equations  of  Motion  for  Submarine  Simulation"  (see  ref  8)  that  would  apply 
to  an  unpropelled  vehicle, 

In  general,  the  terms  which  are  retained  herein  but  are  not  re- 
tained by  either  NSRDC  or  HET/IL  are  done  so  with  the  idea  that  not  all 
submersibles  possess  the  near  fore  and  aft  and  symmetries  that  modern 
military  submarines  and  the  DS.RV  possess. 

Tho  proceodure  to  bo  followed  will  be  to  look  at  the  terms  th  t 
remain  after  the  :  ,■    I  y  torms  have  boon  eliminat*  1  and  tl   1         in- 
cd  by  other  authors  have  been  set  aside, 


3'i- 


Prior  to  lo<  king  at  t!  1 

the  eff<  .  '  ":'  the  choice  of  i     3   poini  f<      '.  , 
invc  '   bed. 

The  most  common  operating  poj    bout  wl  ich  hydrodynamic  fds 

'      expanded  is  some  finite  forward  velocity,  W]  >n  this  is 
done,  the  linear  terms  given  by  equations  (2«5  -  ^3  through  !:r'J   e:  ist. 
These  terms  3      "'  le  the  effect  of  circulation  which  does  not  appear  'in  the 
potential  theory  and  the  Munk  moments  which  arise  from  the  potential  the  or;  -. 
Additionally  these  terms  account  for  the  effect  th  t  some  f:i>  '  '   :'  dtlal 
velocity  has  upon  the  drag  xerms. 

For  the  present  study  of  emergency  ascent  trajectories,  the  ini- 
tial operating  condition  is  for  the  vehicle  to  be  at  rest  in  the  fluid. 
Since  the  force  caui  Lni  th  vehii  "'   to  ascend  is  the  result  of  :■  '' 
ballast  from  any  position  on  the  vehicle,  it  is  as  likely  for  the  ve- 
to start  moving  astern  as  it  is  ahead.  With  this  sort  of  op   ting  cond:: 
tion,  the  most  reasonable  expansion  point  for  the  Taylor  i      '   :'   t)    tl 
zero  velocity  condition.  This,  however,  presents  the  problem  of  complei 
eliminating  all  the  linear  terns  from  the  series  expansion,  since  they  mi  I 
be  evaluated  at  the  expansion  point,  zero  velocity.  For  ex; 


•ft]         ■  ° 

^-    — '  u  =  u. 


since  uQ  -  0 


o 
DispH      the    eliminatio]    of  the  linear  !  3  •  the  Taylor  series  • 

exp    -   3     t,  vq  are  dealin       '  ]      1  fluid  and  the  effect  of  circ     ltion 

will    .  bill    ■     Lse  as  the  vehicle    c         nces  to  move  and  the  potential  theory 

still  indicate      that   t  tk  ]  I         d.st.     T<      ]      Leal  place  fo:     "      -  to 

1      '      Lucl       is,   of  cours    ,    :' n  -;  I  <  '  1  ■  cs  I       ,     '.< 

1 

i   ,  in  the  d   lo]    t  that  "     ,  that  what 
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11;  .'■  ■'-'■"'    i   as  a  line   i  Pfect  is  :\oy   part  of  i        second  order 

off cots. 

p 
The  d]   effects  represented  by  Z  w,  Z  ir J.  etc,  in  the  \ 

^     w         ww 

case  r.iust  now  be  represented  by  only  tl       second  order  terms  'L     w^,    etc, 

WW 

In  the  into:  est  of  d<  .  loping  as  general  a  set  of  equations  as 
possible,  the  linear  terms  usually  appearing  in  the  hydrodynamic  equations 
shall  be  retained,  though  zero  for  this  particular  case. 
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The  axial  force  c  'out  symmetry  terms  is: 

X_  =  X.u  +  X  0a2  -:•  X  r2  +  X  vr  +  X  wq  +  X  pr 
H  .  u    q2     r2     vr     wq     P^ 

+  X  2u2  +  X  2v2  •:-  X  r^2 

u      v      w 

-I-  X.w  +  X.q  +  X  uq  +  X  uw  +  X  vp 
w     q     uq      uw     vp 

+  X  u  v  X  w  +  X  q  (2.5  -  ^9) 

u     w     q 

where  tho  first  two  lines  represent  tho  terms  to  be  arbitrarily  retained. 

Tho  last  two  lines  contain  those  terms  which  require  further  investigation 

before  being  retained  or  rejected , 

The  linear  term  X  ,  vrill  be  dropped  in  favor  of  the  non-linear 

X   which  equally  well  represents  the  drag  phenomena.  This  is  especially 
u 
true  if  consideration  is  taken  of  the  non-dimensionalizing  parameters  in- 

volved.  X  is  non-dimensionalized  by  dividing  by  (■§-^l"U)  where  U  is,  in 

general,  the  velocity  of  the  origin  of • the  body  axes.  Therefore,  rather 

than  leave  the  dimensions  of  X  a  function  of  velocity,  we  can  take  a 

further  derivative  with  respect  to  u  and  eliminate  the  velocity  dependence. 

This  would  then  give  us  the  alternative  non-linear  form  indicated  above. 

The  linear  terms  X  and  X  are  assumed  to  be  zero  on  the  basis 

w     q 

of  experimental  results,   (See  table  I  of  ref  2.) 

The  added  mass  terms  X,,  X,  and  X   are  greatly  dependent  upon 

w   q      uq 

the  vehicle  shape,  If  there  were  fore  and  aft  symmetry  there  would  be  no 

axial. force  resulting  fromX,.  If,  however,  t]    hicle  had  a  t  such 

w 

as  ALVIN,  which  possess  no  fore  and  aft  symmetry,  there  n  3  wel]  "1   forces 

arising  fro  1  X,a  Very  much  similar  arguments  can  be  said         'X 

v  q      uq 

and  therefore  they  h;   1    retained  in  t]         .3. 

3? 


The  second  order  t    X  _  is  assure  3  to  bo  zero  on  ihe  basis  of 

IXVJ 

experimental  results  (see  table  113  of  ref  2). 

Tho  one  remaining  coefficient  X   is  a  second  order  roll  trans- 

vp 

verso  velocity  coupling  coefficient  which  appears  to  bo  essentially  zero. 

II  ~  5«7l.2  }  '  '  Force 

The  lateral  force  equation  without  symmetry  terms  is: 

YtI  =  I.v  •;•  Y.p  •:-  I.r  +  Y  wp  +  Y  wr  +  1    pq  +  Y  qr 
H    v     p1    r     wp     wr      pq      <l? 

-!-  Y  ■/'   •:  Y  vq  +  Y  vw 
y2     vq     vw 

•:•  Y  v  +  Y  p  -!-  Y  r 
v     p^    r 

+  Y  ur  +  I  up  +  Y  uv  (2c5  -  50) 

ur     up    -  uv 

vhore  the  first  three  lines  represent  the  terms  to  be  arbiti'arily  retained. 

On  the  basis  of  the  potential  theory  developed  in  section  II  - 

5 «2,  Y   and  Y   are  of  the  ordor  of  magnitude  of  X.  and  X.  respectively, 
up     ur  v     u 

Strumpf,  in  reference  (2),  retains  these  terms  but  notes  that  there  is 

little  or  no  experimental  data  available  for  Y   while  exper:"     1  results 

for  Y  ^  indicate  that  it  is  important. 

Similarly,  Y   is  also  believed  to  be  of  importance  on  the  basis 
uv 

of  experimental  results  cited  in  reference  (2)„ 

II  -  5.^.3  Normal  Force 

Tho  normal  force  equation  without  symmetry  terms  is: 
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2 

Z     =  Z.w  +  Z.q  +  Z  vp  •:-  Z  rp  +  Z  p 
H    W     q     vp      rp      PP 

+  Z  v2  -I-  Z  v2  •'.-  Z  v2  •:  '/■     v, 
w     w     rr     vr 

+  Z  w  +  Z  q  +  Z  u 

t\T     q     u 

+  Z«u  +  Z  uq  +  Z  wq  +  Z  qq  +  Z  u2  +  Z  ui      (2e5  -  5*) 

u     uq      viq  '    qq      uu      uw 

where  the  first  three  lines  are  the  arbitrarily  retained  terms » 

The  added  mass  terms  Z,,  Z  ,  Z     I  Z   are  retained  on  the 

u   uq   'l  '  qq 

possibility  of  a  less  symmetric  vehicle's  giving  rise  to  this  sort  of  te 
Potential  theory  estimates  Z   to  be  of  order  X,,  which  is  not  negligible, 

Uq  U 

Z,  and  Z   to  be  of  order  X.,  which  was  not  neglected  previously  and  Z 
u      uq  w'  °  J      qq 

of  order  X,  according  to  potential  theory. 

q 

The  Z   term  is  shown  to  be  irapo  bant  on  the  basis  of  experi- 
uu 

mental  results  (see  ref  2). 


II  -  5.^.4  Rollins  Moment 

The  rolling  moment  equation  without  symmetry  terms  is: 

K  =  K,p  +  K.r  +  K.v  +  K  rq  +  K  pq  +  K  vq 
H    p     r     v     rq      pq      vq  ^ 

? 
+  K  \tt  +  K  wv  +  K  wp  -•■  K  _p 
wr     wv     wp  l  21 

-1-  K  v  +  K  p  +  K.  r 
v    p^    r 

-:-  K  uv  +  K  ur  +  K  up  (2.5  -  52) 

uv     ur     up 

where  the  first  three  lines  represent  the  arbitrarily  retained  terras. 

Here  again  there  is  little  or  no  experimental  evidence  i      "    ble 
which  estimates  of  K^,  K   and  K   may  1   i  le.  Pot   bial  t]    ,  in- 
d  sates  that  these  terms  are  small  and,  there:  ore,  ti    1     ,      lectec 
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II  -  5^*5     Pitching  Moment 

The  pitching  .    nt  equation  without  symmetry  terras  :' 

M,T  =.M.q  •:  M.w  +  M  r2  +  M  vr  -:    j   H    p2  •:  M  vp 
II    q    \  rr     vr     rp     PP     vp 

•{-  M  v,2  +  M  wq  +  M  u2 

ww      wq      uu 

+  M  v2  +  M  q2  +  H  ir  ■•  U   q  •:  ]i  u 
vv      qq^     W     q  •    u 

+  M.u  +  M  uw  +  M  uq  (2.5  -  53) 

u     uw      uq  * 

whore  the  first  three  lines  contain  the  arbitrarily  retained  terms, 

Equations  (2,5  -18)  show  that  the  added  mass  tern  M,  is  of  the 

u 

:   3  size  as  X,  amd  M   which  have  been  retained ■  Similarly  M   is  of  the 

q    qw  uq 

same  sice  as  M..  Experimental  evidence  cited  by  Strumpf  supports  the  re- 

w 

tontion  of  M  . 
uq 

The  terra  M   is,  in  accordance- with  equations  (2,5  -  18),  equal 
uw 

to  (X,  -  Z,),  For  most  vehicles  the  Zc  term  is  considerably  larger  than 
u    w  w 

X,  and  therefore  the  term  M   should  ba  retained, 
u  uw 

II  -  5*^6  y-c-.-riria   Moment 

The  yawing  moment  equation  without  symmetry  terras  is: 

N„  =  N.r  •!-  N.p  +  N.v  +  N     pq  +  N     wp  +  M     rq 
tt        r  p  v  pq  wp  rq 

+  M     vnr  -1-  K     vq  -J-  IJ 
wv  vq  vv 

■-   N     r     +  N     •  i 
rr  wr 

-I-  N  v  •!•  N  p  +  N  r 
v  p  r 

+  N    up  •:   N    •.:■  •:  (;  .  ■• 

vp  vu  ru 


l\Q 


where  the  first  four  line;    i       '    Ln  the  terra;    arbii        '  tained. 

Equations  (2„5  ~  IS)   sh  '    1        is  equal  to  I.  »  X., 

vu  v    u   up    vq 

and  N   =  N.,  Since  Y,  is  of  o;  3   of  i     '    ]      ''   t  X,,  N   J   of 
ru    v         v  u   v^ 

order  of  magnitude  large,  i.  n  X,,  N   is  of  order  X.  and  not  ]     ible. 

U    *U  Y 

N.  and  N  i     J  retained  by  the  ot]     thors  a  1  t]  refore,  N   and 
v     vq  up     ru 

should  also  be  ret;  "   1  in  order  to  be  consistent. 


II  ~  6  Equations. ..of  _ Motion, f  or  _  Free_^Ascent 

The  conservation  of  momentum  equations,  equations  (2.3  -  32 
'through  3?) i  £-s  derived  in  section  II  -  3i  the  gravity  forces,  equations 
(2J\-   -  2  through  2,h  -   7)  from  section  II  -  k   and  the  hydrodynamic  force 
equations  developed  in  section  II  »  5  will  now  be  combined  to  give  the 
final  form  for  the  equations  of  motion  for  submersibles  in  six  degrees  of 
freedom  with  varying  mass  and  center  of  mass. 

Several  of  the  terms  retained  are  not  generally  found  in  most 
developments  because  the  terms  are  not  experimentally  or  analytically  ob- 
tained by  present  methods.  » 

As  it  is  not  the  purpose  of  this  paper  to. evaluate  additional 
s  for  the  equations,  these  additional  terms  will  be  retained  but  set 
to  zero  in  the  computation  of  ascent  trajectories. 

The  equations  of  motion  for  a  vehicle  with  varying  mass  will  1 
presented  in  the  following  manner:  the  left  hand  side  of  the  equation 
represents  the  rigid  bedy  dynamics,  the  right  hand  side  represents  the 
hydrodynamic  forces  and  moments  acting  on  the  bod;  and  <  using  the  motions. 

The  hydrodynamic  terms  presented  he   1    ',   n  non-diraensionaliz- 
(  \   in  the  usual  manner  as  described  in  3     nee  (6).  Typical  non  13 
:  "   ^  '-' '   rms  of  the  hyd         i     are  presented  in  the  nomenclature. 


For  :      '  ty,  the  pr:   .  1    1    o  Ltted  :<    '   terras  in  tl   i    '  Jons 
]   ;ented. 

The  equations  of  notion  for  a  freely  ascej  U   .  vehicle  with  vi 
le  mass  are.*' 

AXIAL  FORCE 

in  I  u  -  rv  +  qw  -  x  (q  +  r  )  +  yQ  (pq  -  r)   +  '/Q  Cpr  +  q)J  = 

4  i2  (x  u2  +  x  v2  +  x  v2; 

2     uu     vv     ww 

/°  3 
+  ~  1  (X.u  +  X,w  +  X  uq  +  X  vr  +  X  va) 
£  u     V7     uq      vr     wq 

+  ~  1  (X  q2  +  X  r2  +  X  rp  +  X.q) 
2      qq      rr      rp      o/ 

-  (v)  -  b)  sin  e  .  (2.6  -  i; 

]  \TSRAL  FORCE 

m  |v  -  pw  +  ru  -   y  (r  H  p  j  +  ^-.  (qr  -  p)  +  x^  (qp  +  f  jj  = 

r  1"  (Y  v^  -i-  Y  uv  +  Y  uv  +  Y  w) 

2        VV        V        UV        VW 

P  3 
+  «  r  (Y.v  +  J    ur  -:  Y  ur  +  I  vq) 
6      v     ur      r      vq  x 

P    3 
-J-  o  1  U  up  +  Y  up  +  Y  wp  +  Y  vr) 
^     up     P     wp     wr 

P    h 
+  :l  (Y.p  +  Y.r  +  Y  pq  -I-  Y  qr) 
2      p^    r     pq^1    qr1 

•1-  (W  -  B)  sin  <?>   cor,  0  (2.6  -  2") 


i) 


NORMAL  FOJ 

m  [w  -  qu  •:  p  •  zQ  (p  +  q  )  +  x  (rp  -  q)  +  ;■   (rq  ■:  / ',  |  > 

£  I2  (Z  uu  -:-  Z  u2  +  Z  v2  •:-  Z  w2  •:  Z  v.  ■:  Z  uw) 

2      u      uu      w      -i        I       uw 

+  |  I5  (Z.u  +  Z.w  +  Z  vp  +  Z  vr  +  Z  uq  +  Z  vq) 
Z  u     v     vp      vr      q      wq 

+  §  1   (Z.q  +  Z  p^  +  Z  q^  +  Z  r^  +  Z   rp) 
2      q     PP      qq      rr      rp 

+  (W  -  13)  cos  G  cos  $  (2,6  -  3) 


ROLLING  MOMENT 

I  p  +  I   (r  +  pq)  +  (I   -  I  )  qr 

>        XZ  ZZ  yy   l 

+  m  |yG  (w  +  pv  -  qu)  -  z     (v  +  ru  -  pw)j  = 
~  I3  (K  v2  +  K  uv  •:-  K  vw) 

2      w      v      VV7 

P    h 
•I-  ~  1  (K.v  +  K  up  -!■  K  ur  +  K  vq  +  K  vrp  +  K  vr) 
2     v     p     r     vq     vp     vr 

+  f  l5  (K.f  +  K.*  +  Kppp2  +  ?pqpq  +  Kqrqr) 

+  (y  W  -  y  B)  cos  6  cos  $  -   (z.  W  -  &JB)  cos  6  sin  d 


(2.6  -  /-!) 


1-3 


PITCID.J   '  1 

I     4  +  1       (r2  -:•  v)  +  (]        -  I     )  rp 

yyH  xz  XX  zz 

+  m  fzG    (-a  +  qw  -  rv)  -  x     (w  +  pv  -  qu)J   - 

~  "3?  (M     u2  +  M  u2  •:   M     v2  +  M    w2  +  M     uu  +  j    uw) 
2  uu  u  w  uw  i 

+  ~  3      (M.u  +  M.w  +  M     vp  +  W     vr  +  M     «q) 
2  u  w  vp  vr  wq 

+  «  1     (M  uq  +  H     vq) 

^  q  uq 

+  1 15  (,;,q  +  lippP2  +  i-iqqq2  +  Mr/2  +.MrprP) 

-  (xqW  -  XgB)  cos  0  cos  $>  -  (zG"W  =  z_B)  sin  6  (2.6  -  5) 

IAWING  MOMENT 

1  r  +  I   (p  +  rq)  +  (I   -  I   )  pq 

zz    xz  x  yy       xx    * x 

+  n  Jx  (v  •!  ru  -  pw)  -  yu  (i*  +  qw  -  rv)  = 

~  1^  (N  v2  +  H  uv  +  N  uv  +  N  w 
2      w      v      uv      vw 

+  V  1*  (NiV  +  N  wr  +  N  wp  +  N  ur  +  N  ur) 
2     v    wr     vp     i''     ur 

>°  4  / 

+  -  1  (W  up  +  N  up  +  N  vq) 
2      up      p      vq  1 

+  r  1?   (N.p  +  N.r  +  N  r  +  M  pq  +  N  qr 
2      px    r     rr      pqr^    qr1 

+  (xjjW  -XgB)  cos  0  sin  0  +  (yj  -  yB)  sin  0     (2.6  -  6) 


/;': 


13  •-  7  S    ry; 

i    Li  syi  bems  to  bo  used  in  the  ]   '        discussed  and  t] 
transfoi   bions  from  the  earth  fixed  to  the  body  fixed  axes  s  'e  d  veloped. 
A  derivation  of  th<  d;    '   "'  equations  ;   a  vehicle  with  a  v<  3 
ss  and  center  of  mass  is  then  made  using  a  Lagrangian  formalism. 
assumption  that  the  reduction  in  mass  due  to  releasing  ballast  can  be  re] 
sented  as  a  quasi-steady  process  is  discussed* 

The  forces  acting  on  the  body  are  then  developed  by  expanding  a 
functional  repr(  '  lion  of  the  forces  in  a  Taylor  series.  The  terms  in 
the  series  arc  then  discussed. 

The  equations  are  then  non-dimensionalized  for  final  presentation. 
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CHAPTER  III 
SOLUTION  OF  TKE  EQUATIONS  OF  MOT 

In  the  previous  chapter  a  system  of  equations  of  motion  were  de- 
veloped for  which  a  solution  must  be  found  in  order  to  obtain  the  ascent 
trajectory  of  a  vehicle.  It  is  possible  that  an  analytic  solution  to  this 
set  of  simultaneous  non- linear  difJ     I.ial  equations  could  be  found,  he 
ever,  it  is  much  more  practical  to  assume  that  a  stepwise  linear  approx." 
tion  to  these  equations.  'This  approach  has  been  used  in  many  simulat' 
studies  of  which  the  DSRV  control  system  simuD      '  see  ref   5)  and  deter- 
mining the  effect  of  hull  shape  non-li.     ies  (see  ref  4)  are  jr  b  two. 

The  equations  .      .loped  in  C]  }  I   II  conform  to  the  notation 
used  by  Strumpf  (see  ref  2)  vrhich  does  not  conform  with  either  the  hydro- 
dynamic  terms  format  used  by  the  HIT  Instrumentation  Laboratory  for  the 
DSRV  control  system  studies  and  simulation  (;     '  5)»  or  the  t  ! 
equations  of  mo'/    For  submarine  simulation  as  used  by  NSRDC  (see  ref  8), 
In  order  that  a-'  L]  ble  hydrodynamic  coefficients  be  utilised,  the  notation 
of  Chapter  II  will  be  modified  in  this  chapter  to  conform  with  that  of  the 
NSRDC  equations.  The  choice  of  the  NSRDC  form  over  the  KJT/IL  form  was 
made  for  the  reason  that  it  is  more  likely  that  vehicle  coefficients  will 
be  obtained  from  KSRD3  than  from  KIT/IL  which  must  also  get  its  data  fi  - 
other  sources.  The  KIT/IL  notation  does  possess  the  advantage  of  ha\ ' 
the  e'   tsional  forms  of  the  coefficients  independent  of  velocity. 

Before  a  stepwise  solution  to  the  equations  i  hied,  the 

be  put  into  .  I  .•     bio  to  t      hnic  le. 


*i6 


-1  '  '  "*  1  fiovii  sod,  Equations ip  of ,  jjotj-On 

Tho  vehicle  equation!  i   motion  are  generally  c  the 


form: 


mI!=-n(l  *»+Wx(v+1  :  KH+t^+f^  (3.1  -D 


I:  ( 


W  x  IW  -  mBG  x  (^  +  W  x  VJ  •;  M^  -J-  HfiFF     (3.1  -  2) 


dV  _  •   • 


d  . 


where  ~-  =  u,  v,  w.  —  =  p.  q,  r  and  I  is  tho  inertial  tensor. 
dt  dt 

To  determine  ascent  trajectories  in  a  stepwise  linear  fashion 
on  a  digital  computer  it  is  useful  to  rearrange  the  above  equations  so  that 
all  the  derivatives  are  on  the  left  side  of  the  equations.  The  resulting 
equations  are  of  the  form: 


dV' 
dt 

d£' 
dt 


*x<v  +  wx»)  +  r   +fot 


(3.1  -  3) 


— »-     '.:_t. 


vhero 


W  x  !••.  •  mBG  x  (W  x  V)  +  BG  x  Wy  +  Hfm)  +  K^.  (3.1  -  fc) 


dt 

.   Ftt"! 

St 

dt 

~   lMJ 

-Ok 

dW 
dt 

(3.1  -  5) 


and 


|i*j  is  the  six  by  six  derivative  coefficient  matrix  given  on  page  /,'8. 
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dV 
The  actual  deri  -  and  •       ;        obi  d  frora: 

dt 


■ 


( 
dt 


[  ]  ' 


'<.-. 


(3.1  -  6) 


whe 


[m]-i  5 


is  the  inverse  of 


H 


clV    „,~^ 


Tho  equations  to  be  solved  j   '  '  and  -    £ ,  'e  then  t  '■    by  the 

dt  dt 


f OIIOWJ        ! 


dt 


=  -  ra   j«  rv  -!•  c.    •    :   ,    (q     •:    :>/-)  •',-  yGPq  +  sIqP^I 


+  C;   I2    (X      U2    +  X      V2  +  X      v?) 

c  uu 


V  •, 


\  . 


P  -,; 


+    ~    \J     (X        Uq    +  .X        VI-    +    X        1  q) 


uq 


\. 


\  \ 


•:  • 


P  J 


j  ''   (X    o2  +  X     r2  +  X     rp) 
qq  x  ) ,  rp  ^ 


(..   •    B)  sin  6  +  T 

bFF 


d]  ' 
d1 


Dl 


[-  pw  •:    ru  -  yQ   (r2  +  p2)  +  zGqr  +  xQqpl 


(3.1  -  7) 


+  ~  l2  (I  ,   ,v  I  (v2  +  w2)  * 1  +  y  uv  +  y    w  •:-  i,u2) 

<c  v|v|      '  '  V  vw  * 

+  I  I3  (X  up  +  y  ur  +  Y     vq  +  Y     vp  +  Y     wr) 
c,            p            r  ""^^  ■*r'' 


£  n3  , 
2  J- 

r  ,4 


vq  wp 


Wr 


(Y        -yvr     (v'  +  w*)  2  |r|) 
v|r|  |v|  '     « 


+  l6  1"    (Y     pq  -:-  Y     qr)  +  (W  -  D)   sin  e*  cos  0  •;    3 

_  .        Pq  ci:i-' 

(3.1 


8) 


' 


d  ' 


*^2(Z,     u|w|+Z       w|(V 

2  VI 


1  . 


cv.< 


+ ^  <V>  +  V"  V1 :  V^i  |(v* ':  A  ' llql) 

*  i  3'!' ( v' +  v2  +  v2  +  Vp) 


+  (W  -  B)  cos  0  cos  j$  +  ZFR, 


. 


(3.1  -  9) 


I     )  qr  -  m  ["-  z     (-  pw  +  rq)  +  yQ  (pv  -  qu)J 


-  I     pq  +  (I 

+  £  I3  (K   ,    ,v  |(v2  +  w2)  tj  +  Kuv  +  K     vu  +  K»u  ) 
2  v|v|  * 


+  ^l5  (K    plpj  +  K  pq  +  K  qr) 
2      p|pj  I  *    pq      qr 


+  Cy  W  -  yBB)  cos  0  cos  $   -  (zGW  -  2£B)  cos  0  sin  $   +  K^ 

(3.1  -  10) 


:    -  I       (r2  -  P2)  +  (I       »I     )  rp 

XZ  ^2 


I 


('    G    C-qu  v  p   )  H    Zg    (q     ••  rv)J 


^3/„     2 
i  v 
u 


2: 

v     +  M  I     ,i 
w  v  w 


+  §  3/  (H_u&  +  M_    vc  +  M        ,i     |(v     + 


v 


/?n3    /„ 


+  ~  I"*  (  .    (v"  +  \r)  ' 

£  V 


I   M.    ,u|v|) 


1    (       vp  •:  M    vr  •:  M  uq  •:  M.   ,     Kv2  •:  w2)  2  q) 
vp  vr  q  ^  |w|q  '  '  * 

r  1J  (M     p2  +  M    ,    ,q|q|  +  M     r2  +  M     rp) 
2  PP  q|q|     I*'         rr  rp 

-  (xJW  -  XdB)  cos  G  cos  $  -  ("a-W  -  z  B)   sin  0  +  2-1 .,  ,_ 

u  1>  u  c  i^fr 


i°5 


(3a  -  n) 


dt 


I     rq  +  (1       -  I     )  pq  -  r.i     -  yp   (-rv  +  qv)  +  xu   (ru  -  pw) 

zx  xx         yy        "  L-      C-  u  J 


P    3 


yy 

.2    i   __2  \   2 


■   «  ..2 


{-  ~  l-'  (N        ,v     (Vs  +  w  )  2     +  N     vw  +  Nuv  +  l\Ur) 
2  vjv|  '         vw  V 


+  £l     (N     vr  +  N     vp  +  1"     vq) 
2  wr  vp  vq 


+  "J  1     (N  ur  +  N  .    . 
2  r  ■ 


(v2  v  Vr2)  *  I  r  +  N  up) 


+  ~  l5   (M    ,    .r|r  |+  II     pq  +  N     qr) 
2  r|r|  pq^         qr* 

+  (x  W  -  XgB)  cos  0  sin  <f>  -V   (yQW  -  y  B)   sin  6  +  K 


(3.1   -  12) 


In  these  equations,  unlike  equations  (2,6  -  1   through  ?)  from 

which  they  come,   the  linear  terms  I  ,  Y  ,  Y  ,  Z   ,  Z   .  ]     ,  ]     .        ,  ]     .   N 

v   p   r'  u   w   u  vr      q   v   p 

and  N  are  combined  vith  the  non- linear  terms  Y  ,  Y  ,  Y  ,  Z  .  Z  ,  ] 

r  uv   up'  ur   uu   uw1  uu 

J-,  .»  J:^p  N  i  N  ,  N  ,  This  is  done  for  the  reason  that  1 
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inseperab]  frc  i  t    i  Dde3  test  data  reduction  \     int.  'j   I    :  Z  ., 

r  l 

Z  ,  M    and  i'J   he-.  1         ■  o  as  to  include  all  of  the  KSRDC  ^   ffi- 
ww   I  w  I     * 

c-  nts.  The  ten  s  X^,  X^,,  Z£FF,  K^,  M^,  and  ^   ]    1   n  included 
in  order  to  allow  for  the  possibility  of  there  being  some  small  effector 
forces  on  the  vehicle.  This  inclusion  in  no  \  ,  affects  the  '.'   gitimacy  of 
the  equations  so  long  as  the  force  does  not  cause  the  pre  h    '  ance  of  one 
velocity  component  . 

With  the  equations  in  this  form,  a  stepwise  linear  solution  can 
be  developed  which  will  be  capable  of  being  programmed  for  use  on  a  high 
speed  digital  computer. 

Ill  -  2  Stepwise  linear  Solution 

A  step. rise  linear  solution  is  one  in  which  the  accelerations  are 
taken  as  constant  over  a  given  time  interval.  The  accelerations  that  exist 
over  a  time  interval  are  determined  from  the  velocities  that  existed  at  the 
end  of  the  previous  time  interval  and  the  total  weight  removed  from  the  ve~ 
hide.  Thus,  the  right  hand  side  of  equations  (3.1  -  3  and  3*1  -  ll)   can  bo 
determined  at  the  beginning  of  each  time  step.  The  inertia  ten   '   : 
must  also  be  recomputed  before  each  step  since  th     s  is  changing  in  a 
stepwise  linear  fashion  as  described  in  section  II  -  J,     With  the  foregoing 
information  equation  (3.1  -  6)  can  bo  solved  for  the  body  axes  acceleratio 

The  velocity  of  the  vehicle  in  the  body  axis  sysi/    y  then  bo 
found  by: 


V 


W 


dt 


(3.2  -  1) 
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wh  'c  t  5      '    at  which  the  ascent  sta;  '       b  is  the  present  t: 
o 

The  inl    1  ma;   1   over,  be  repres     I  <    s  ■    each 

step  i  bepi-d.se  solution,  Thus: 


v 


V 


n  =  1  t 


Dj 


dt 


J  n  -  1 


(3.2  -  2) 


v1  re  t     ;  d  t  are  respectively  the  times  at  the  begin:  '   and  end  of 
n  -  1      Y> 

the  interval.  The  solution  I  •  been  specified  as  si   .  se  linear  and  tj 

accelerations  are  to  be  obtained  as  constants  for  the  duration  of  an  int 

val,  therefore,  the  accelerations  can  be  removed  from  the  ini  gral.  leavii  g 

only  the  trivial  integration  of  dt  from  t    ,  to  t  0  Defini]  g  a1  as 
"  °  n  -  1    n 

t-  -  t    „  ,  equation  (3*2  -  2)  becomes: 
n    n  -"  1 


V 


W 


Jt 


n  =  i 


i- 


At 


(3.2  -  V 


cn  -  i 

which  gives  the  translational  and  rotational  velocities  of  the  vehiclo  in 
the  body  axis  syst< 

In  order  to  obtain  the  actual  position  of  the  vehicle  relative 
to  its  inertial  starting  point  the  velocities  obtained  above  must  bo  trans- 
formed into  the  inertial  axis  system  by  uue  of  the  inverse  of  i       1    for- 
mation matrixes  developed  in  section  II  -  2,3«  The  vole  '     Ic,  y,  z,  $, 

•   * 

G,  y  are  obtained  from: 


=  T. 


J. 


u 

8 

W 

' 


•  1 


p 


(3.2  -  ;:); 
(3.2  •  5) 


The  position  and  or:      on  in  the  inertial  axis  system  is  then 
detr.     d  by  integrating  the  respective  velocities  from  t  to  t.  Here  a- 
gain  the  solution  is  stepwise  and  may  bo  represented  as  a  sum  of  the  indi- 
vidual steps.  In  this  case,  however,  the  integrand  x,-,,  yu,,  etc.  is  not  a 
constant.  It  is,  instead,  a  linear  function  of  time  since  the  acceleration 
from  whence  it  came  is  a  constant.  Therefore  the  <.  j  '  Ion  describing  the 
vehicle  trajectory  and  orientation  becomes: 


>■ 


: 
E 

h  ( 

~i 

_ 

?■-• 


XE 


r 


;•• 


(3«2  -  6) 


>■ 


\jhe 


re  Xg  =  (xg,  yE,  z£)  and  ^  =  (^  C,f) 


The  computation  of  x_,  and  (   are  tl  i  final  computations  of  the 
step.  The  procedure  is  then  repeated  until  some  maximum  time  is  read 
or  some  predefined  position  is  reached. 

In  oidcr  to  verify  the  result  obtained  by  the  foregoing  pro- 
cedure, and  independent  trajectory  determination  must  bo  made. 

-'--I--,-  "  3  Pfls^imensional,  Ascent.. Tr; 

Now  that  a  three  dimensional  second  i  "   i  thod  of  determining 
free  ascent  trajectories  has  been  obtained,  some  check  on  the  results  of 
this  procedure  is  in  order.  The  most  1  gical         :'       duce  tl 
problem  to  one  in  tho    I  Leal  plane  only  and  thereby  ,  I         com- 
plexit;  of  tho  prob]     If  a  c   '   lion  of  the  equ 
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made  at  t!        the  resultii       Lons  would  bo  no  different  from  t: 
previously  obtained  except  that  some  of  the  cooffic:'  ro. 

In  addition,  the  solution  of  the  two  dimensio  1  eq 
trajectory  would  rcquiro  the  aid  of  a  high  speed  digital  computer. 

To  avoid  the  necessity  of  employing  a  compul   ,  •  :  implified 
method  described  by  Giddings  and  Louis  in  reference  (9)  will  be  used.  This 
method,  as  applied  to  the  problem  under  consideration,  reduces  to  a  one  di- 
mensional solution  of  the  equations  of  motion  for  a  vehicle  with  varying 
mass. 

The  one  dimensional  solution  determines  the  vertical  position, 
velocity  and  acceleration  of  the  vehicle  during  its  ascent.  This  informa- 
tion provides  a  sufficient  check  of  the  results  of  the  computer  solution  of 
the  equations  of  motion  developed  in  Chapter  II.  A  description  of  the  one 
dimensional  solution  follows. 

Defining : 

B  -  the  net  buoyant  force 
b  -  the  time  rate  of  change  of  B 

F  -  the  sun  of  all  forces  in  the  vertical  direction 
g  -  the  acceleration  due  to  gravity 
k  -  the  virtual  mass  coefficient 

m*  -  the  mass  of  the  vehicle  plus  the  virtual  mass 
W  -  tl   :'".'   i1   oous  vehicle  weight. 

z,   is,  z-  the  vertical  displacement,  velocity  and  acceleration  respec- 
tively 

Z  t  Z,  -  the  croi   low  d:   ;  and  ;       ss  coeffi'  Lents  :     ' 
ww   w 

Nexrtons  laws  of  motion  can  be  i       I  as: 
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F  =  m'z  (3.3  -  1) 

The  force  F  is  made  up  of  the  w«  '  ;]  '  of  the  jettison        st, 
B,  plus  or  minus  the  hull  drag,  B,  depending  on  the  direction  of  motii 
The  force  due  to  the  jettisonning  of  ballast  is  given  by : 

t 
B  =Vb  (t)  dt 

o  (3*3  -  2) 

Letting  the  integral  be  represented  by  a  series  of  finite  steps  the  buoy- 
ant force  becomes: 


N 


B  ^   >   b 
n 

n^l  (3.3  -  3) 

where  the  b  are  a  sequence  of  finite  weights  to  bo  jettison     id  n  is 
n 

the  number  of  intervals  ellapsed  since  time  zero. 

The  force  due  to  hull  drag  is  represented  by: 

B  -§/>l2  Z  z|z|  (3.3  -  k) 

WW 

The  mass,  3a',  of  the  vehicle  can  be  written  c 

W 
m*  =  -  k  (3„3  -  5) 

where    k  =  [z;  (£/> 13;  /  (W/g)J  +  1  (3.3  -  6) 

Substituting  these  exprc      '       :  into  '       1  (3*3  -  i)  the 

tion  oi        bion  bee 


(Z,  (i^l3)  +  -)'  -  Z   (.  p  I2)  'M\   .   }~   b  =0   (3,3  -  V) 

Tho  technique  of "  step  by  stop  inl     Lon  sir/''    Lo  that  used  in  section 
III  -  2  may  now  be  applied  to  this  equation  in  ore  one  di- . 

nencional  trajectory, 

III  -  h     Sum'.'ary 

The  notatis   and  terms  format  to  be  us€     the  computer  simula- 
tion is  discussed e  The  notation  used  in  the     bions  of  :  tion  is  modi- 
fied to  conform  with  that  of  KSRDC  and  the  equat.'    ;   re  rranged  for  d 
ital  computer  solution. 

A  solution  method  is  developed  utilizing  t    stej  Lse  linear  tech- 
nique o  It  accepts  the  vehicle  velocities,  position,  o:         and  buoy- 
ancy as  initial  conditions  and  using  equations  (3*1  -  6  through  3ei  -  12) 
it  computes  the  accelerations,  velocities,  and  disp*'      bi  oJ  bhe  vehicle 
after  a  time  interval  At.  Before  each  stop  the  weight,  CG,  mass  and  mora 
of  inertia  are  adjust .  3 

Finally  a  simplified  one  dimensiona]  i  '."     3  is  devised  for  compar- 
ison with  the  computer  results. 


5? 


( 

In  Chapter  II  the     tions  of  motion  were  derived  for  t 
ascending  vehicle  with        .■  iss  and  center  of  i    .-  A  ste] 
solution  suit  ":    for  programming  on  a  digital  cc  d  eloped  in 

Chapter  III,  and  the  actual  progr<   is  presented  in  Apj  ndix  A.  An  equation 
for  a  one  dimensional  ascent  trajectory  was  also  developed  in  Chapter  II  to 
provide  a  check  for  the  three  d*Jb  ansional,  six  degree  of  freedom,  ascent 
trajectories  pro;  ..■■  ,  A  program  to  solve  this  i    1:  tial  «  .*  a- 

tion  is  presented  in  Appendix  D, 

In  this  chpater,  .the  computer  simulations  and  results  of  these 
simulations  are  discussed,  conclusions  drawn  and  '    for  future 

work  in  this  area  made. 

IV  -  1  Re  suit  s_of  _C  omput  er  ,s  ii"1 1 "'    ns 

Computer  simulations  of  the  ascent  t]  j   I   ies  for  i  .  1  -le 
lar  to  the  Deep  Submergence  Rescue  Vehicle  were  conducted  v  '  ,  the 
IB  36O  co  puter  of  the  MET  Information  Processing  Center.  These  simula- 
te is  were  conducted  pr:    "  \.  ;    the  purpose  of  del  ,  Lng  the  ascent  tra- 
jectories program  and  determining  its  oj>erati.      racteristics  and  second- 
arily to  study  the  motion  tendancies  of  the  DSRV. 

IV  -  1.  ,1   Prog:    '     '         '   Lstics 

The  ascc  '  I    ajecto]  '   ]           Iculating  tl  '  Lcle  trajec- 

torii   -1  ]   3  upon  the  ass    I   time  incre   t  f  or  ii  y.  In 

V      ['        '        j     1  can  be  cl    \            '  '          '        3  ' 

all  right,  Noi   that  1   latl   •  '    .     "all  t  ", 


of  just  t\       ■ ."  ;]  I  '•      it  requires  a  c'     Lon    ;      tors 
influ   "    tl   siti  " '■  n  and,  therefore,  the  former  c      '   •  111  be 
d«  '  :  "   ji  fi3  :  I  . 

Should  V       '<      3  increment  cho    1   too  ':    ,  acting 

Oil  th   i   Lcle  will  not  be  damped  in  s   \  L  fash:  i.  A  vehicle  und 

the  infli     of  j  si   3y  force  vn.ll  contir   to        '   unti]  this 
force  is  overco   b;     ''    force,  and  tl     puter  does  not  see  another 
force  until  I  he  p:      '..'■   '       t  ends  and  the  forces  result:   from 
tl  3  present  motion  can  be  calculated.  If  the  time  increment  is  too  long, 
the  motions  calculated  by  the  computer  bee     cessive.  The  eventual  re- 
sult is  that  the  cc  ]  ter  run  is  3    turely  terminated.  "■  1     tual  cause 
of  terrain  tion  (  )  be  either  pri  ;r,  1  <.  .      ed  o:  i   hine  caused. 

A  run  is  terminated  .by  th  pro^ tb  time  limit  or 

tl   depth  ";  "  ''  is  exceeded,  Since  ]    '     tl   I'   Id  'it  or  re; 
depth  limit  in  a  natural  fashion  do  not  constitute  being  ]    - 'are,  they 
are  not  consider     re,  The  deptl  lii  Lt  c   ,  however,  be  :•   hed  in  an 
unnatural  fashion  as  is  indicated  in  Appendix  C,  section  C  -  3< 

A  run  is  terminated  by  the  com]  iter  if  the  n  1    of  exi: 
vali   calculi  ' "'  ns  becomes  excessive, 

Te:>  "   tion  for  the  latter  cause  occurs  before  '       1  outpul 
1   '  ■  3  and,  therefore,  the  only  printed  output  recei^  I  wil3  1 
init:  ''   <  ■     "tion  printout.  Termination  fc  tfo  f     cause  : 
■•  ]  '■    no   1  1     '  pui    in  printed . 

lv hen  the  time  inc]      .     d  is  1        ; 
off  e:    d     '   th  actual  calculations    it!   output  2  '      less. 

Th:  i  then  ]         ' 
'  t   I  !  :  '      t" .  Fact*    hicl      '        '    ted  bj 


clue  val  rate,  vehicu]  ■•  velocities  and  ;    .'       , 

Lular  stability  both  static 

The  ballast  ro"'     :     '  (   1  time)  can  be 

assumed  to  be  fixed  for  i    given  vehicle  since,  in  emergency  conditions,  ■■' 
ballast  will,  in  g«   ,  1,  be  role,         high  a  rat  i   ]  .  ;ible.  Even 
if  this  is  not  done  the  ballast  release   ;   vould  normally  be  specified  and 
not  left  to  the  needs  of  the  simulation .  Further,  the  program  in  no  \ 
affects  the  ballast  release  rate. 

The  velocities  and  accelerations,  o:  t]   other  ha:  \   are  directly 
affected  by  the  time  increment  used.  The  change  in  velocity  is  determined 
by  as j   Lng  that  the  acceleration  is  constant  over  the  time  period  and  is, 
therefore,  <  ;.  1  to  the  product  of  the  time  increi   I  and  the  acceleration. 
The  i    ccle-j  "'•ion  existing  during  a  time  increment  is  dependent  upon  the 
forces  acting  on  the  body  during  that  increment.  The  hydrodynamic  forces 
are,  in  turn,  dependent  upon  the  velocities  computed  di  '    Lb    svious 
increment.  This  interd  ;   lance  of  forces,  accelerations  and  velccit:' 
can  be  the  cause  of  p:    bure  to:*  '   bion  if  th  '. '   inert    '  .' 

ge.  For  instance,  a  3-    inc:    nt  vould  cause  the  acceleration  to  act 
for  too  long  c    Lng  tl   veli  it;  bo  1        "   ,  \ '  ,h  in  turn  would  c 
the  forces  in  the  next  increment  to  be  excessive,  etc.... 

The  statj   stal  Llity  of  tl   srehicl   '   dependent  only  upon  t] 
shift  of  the  CG  :  '  bive  to  tl   C  ,  :    \   :'  ■  ,  '    F<   ,  a  funct: 
ballast  re]    \   rate  only.  T]   d;    i  :  '  ■  1  '  :  3  ty  is,  h  ,  de] 

i  the  vehicular  velocities  and  ac<  "'   tions  i   well  i      V.      cl     :' 
i    .   ;ed  by  tl   reduction  i 

The        of  this  c      '    Ls  that,  t]    '    "      it  to  1 
used  fo]   vehicle  j  will    dictated 


1  ise  rai     It)    pected     '    tions.  Th 

•  ill  indie  be  t]   order  of  i    Ltuc    11      '    '      ."'  ' 
bl   :'       -b. 

With  ;"   proper  tj   3  ic:    nts  cho    I   j   '   ies  can  be  co. 
puled  for  a  variety  of  initial  cc;      is, 

IV  -  1 .2  DSRV  Motion  Tend; 

The  ascent  trajectories  generated  b;  the  i  ]  not 

those  th  b   uld  be  followed  by  the  DSRV,  since  the  coefficients  used  in 
the  progra  -  come  frc  i  DSRV  model  tests  with  the  propeller  running.  They  do, 
however,  provide  an  indication  of  the  ascent  trajectories  and  motion  ten- 
dancies  of  that  vehicle.  Simulations  were  made  with  the  body  ,  :*  dtially 

coi  '  '  :  '.  •  :'  bh  the  inertial  axes  and  with  an  in:  '  '  ]  :■  oil    ]   :*  ^osed. 

The  i        Its  of  the  si  "'  '  '  .    ith  the  body  i  :'  '  '  '  •  lly  coin- 

cident with  the  inertial  axes  clearly  indicate  the  coupling  of  pitch  ,■ 
surge  which  bring  about  th  traji  :torie:  oi  .  .'  u   IV,  1. ,  This  motion  I 

s  even  more  cl<  i  in  the  velocity  plot'  of  Fj     31V, 2.  A  tend   :;  to- 
d  t    .:  .!'■   pitch  angle  when  dumping  both  tr"   '•  nks  into  the  ]     .oir 
is  also  demonst  -  '  1, 

The  simulations  involvi  g  ;   initial  roll  angle  i]  "'   '   bhe 
s  ill  amoi  n1  of  roll  damp:   -    '     ■"'""'  :.      Lated  w: ve- 
hicle, They  also  3  '   be  a  slight  :  '    force  roll  co   ''      is  evidenc- 
ed in  tl     '  '  ■  t  in  Appendix  C,  i      '  '■  n-  ■■  2,     T3   ]  k  of  stal 

Lf  by  c    Li   bl         bo  I   :   '       '   ly  wl 
1  '  ;  a  time  increment  tl    ,     I  si      'ul  for  the  zero  rol  I 

b  of     tenth  that  i     Lthout 
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IV  -  2  C  onclusions  an '  Lions 

'  :  Lmulat:      '    \  ■  terisl  cs  ha\ 

onstrated  the  wort]     this  pro;  r   in  del     '    I     '       a  vel  Lcle 
as  it  ascends  under  the  influence  of  ":   .  ncy  alone.  They  ]     she    thai 
a  vehicle  such  as  the  DSRV  is  qi       sitive  to  roll  without  the  stabiliz- 
ing effect  of  a  relatively  ]     Fo    d  velocity. 

The  v  M      vehicle  du  5    ■  scent  indi 

U  I  tl   e  is  a  predominant  velocity  cc.     it  about  which  ve   could  expa 
a  Taylor  series  once  the  initial  acceleration  ph    •   : '   ■      ]    been 
passed.  The  most  likely  time  for  a  shift  to  this  sort  of  ..  ■         is 
I   nty  seconds  aft  c*  debal]    ■    ..   1    1   n  completed.  At  th,  I  :     '    ve~ 
] I  3     ears  to  h    r<   hed  a  terminal  conditj   in      u 
bch     bei  tial!   ■    I   t  limits 

The  simulations  h  3   ted  that  a  vehic!    '        ]     ct 

:    and  starb   d  symmetry     '   s  no  mo: 
of  motion      tl  re  are  no  sid  -  :  ces  present. 

It  is  :       led  that  future  studies  1     le  using    r'ficiei 
obtai]  if]   i    lei  '.   ts  conducted  without  ]     '    _.  running.  Also,  in 
vie:-  of  the  slow  ascent  velocities,  it  may  pr<      thwhile  to  ]    Li 
equat:    to  include  the  effects  of  an  ascent      lso   hici       ove 
necessary  in  order  to  get  s   more  rapid  as   :-  in  other  th  y  sit- 

Mens. 
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A 
C<       UTER  PROGRAi-1  FOR  ASCENT  TRAJ  IES 

/    •    1     General t 

The  computer  program  to  i    ilate  f         1   '     Les  is 
writ!  •  3^0  -  65  digital  compul  2tj.cn  Proc- 

essing Center.  The  program  is  written  in  (see  ref  11  and  12), 

The         Ices  xise  of  the  matrix  ines  contained  in 

vers:  i   "Scientific  Subroutine  Package"  (ref  1,0  and  the 

math   bical  functions  contained  in  the  standard  "Library  Subroutines"  ( 
14). 

The  system  subrouti  : ; 

•rary  Subroutines ; 

C  S  -  computes  the  cosine  of  an  angle 

iputes  i '   i  mgent  of  an 

SQl'iT  -  c(    I  t]   square  re  '    •    ber 

ABS  ~  cc  ]  '  5  the  absolve'    "'   -of  a  n 
Scientific  Sub.   '     ; 

!  PRD  ~  comj  '  :  the  ]   '  ct  of  t  0  m   '    ' 

I'UNV  -compul  the  :       of  1 

Extensive  use  .'     le  of  the  NAKBLIST  feal    of        IV 
input  and  output. 

The  program  is  parts s  1 

SI  :0UTII\'ES, 


i    -1,1  MAIN  Pro 

The  1  '"'••  i      -  acts  as  an  executive  reading  in  initial  c 
[licle  par,    rs  and  program  contro]  ;  i  ailing  the  var- 

ious subroul        led  to  sol-.   I       '  '  ons  of  motic  '  the 

j        '  the  c      tionsc 

The  data  read  in  is  contained  in  fj       LISTS,  1  ,   FLUID, 

,  '■        L  (     ,      le  controlled  f<    b  '  .    t- 

put  is  again  in  ace   "  ce  with  the  N.        plus  a  s: 
POSIT,  co]  baining  "  .  ■  '  Lod. 

Part  of  the  output  from        a  set    puched  i  on- 

.1  all  the  p;    bers        SIIST  INCOND,       output 
input  data  for  another  run  in  order  to  continue  the  same  trajectory  if  it 
is  desired   A    Lraum  of  steps 

A  -  1 ,2  Sul    '  '   ' 


COSFF  reads  in  the  non-di]    onal  coej  '  ly 

pri  '   :'    Tor  output.  The  subroutines  then  cc . 
factors,  diraensionalizes  the  coefficients  for  use  in  tl 
and  prints  the  •"'    sionalized  c(  i  i<    '      output. 

The  coefficients  are  read  in  and  pr.       '    I 

ci  ;  ?s. 

The  diraem  '•     I    ss  ccsffii  '  by  1 

dimension;  L 

A  -   1.3     Sub 

EALAST  cc-   '  :  the  changes  in  weight,  mass,  cenl   of 
end  ;    its  of  inertia  due  to  ;  :ts 

sul   >ut:  i  is  c: 


6( 


trajectoj .      debal  place. 

BAIAST  calls  subroul  '•  bed  in 

AHATRX  varies  only  <      bl   del     bing  phase  of  the  ascent. 

A  maximum  of  twenty  3   ' '    for  1  1]   b  weighi 

There  is  no  printed  output  from   I  '  T« 

A  -  1.4  Subroutine  / 

AMATRX  computes  the  six  by  six  ma1  i   ;   acceleration  te:     3f- 
eficients,  AM,  and  then  inverts  the  array  for  use  in  the  solution  of  the 
equations  of  motion. 

A  -  1.5  Subroutine  HYDRO 

FYDRG  computes  the  hydrodynamic  forcing  terms  and  the  g:» 
forcing  terms  and  then  sums  them,  plus  any  ef J        ces,  to  get  the 
total  forcing  terms  for  the  equations  of  motion, 

HYDRO  uses  the  velocities  and  anglei  resulting 
stops  in  order  to  co:  p  ite  the  hydrodynamic  a  \   gravity  force  cc.     '  ■    act- 
ing on  t   vehicle  during  the  present  step* 

There  is  no  printed  output  from  HYDRO. 

A  -  1 .6  Subroutine,.  TIUJEC 

TRAJEC  computes  the  trai  "!  '    ]       '      1  velocities  of  the 
vehicle  in  both  body  fixe  :     earth  :'.'   '  -  '        ;.  In  cider  to  do 
this,  TRAJEC  also  computes  tl   i    !  '  Lonal  and    u]  '     'or- 

I  bion  ]  itrices  '''    ■    and  TINV. 

If  the  pitch      becomes  ninety  deg     the 
trix,  AINV,  vriDJ  blow  up  (■   ]    b:   lly  s;    '  g).  At  '    tine,  the  pro- 
!     '  :    -'  '  a  mess 

■tor  systc  •  •    apt,  IV 


during  corapul      of  the  bi 

cosine  of  the  pitch  i    gle.  At  this  point  the  ..         bake  the  . 

corrective  action  of  assigning  \  Lue  of  10''::'69  to  each  c  nts. 

A  "  1.7  Subroutines  POSn 

computes  the  time  ellapsed  and  the  position  of  th«  i   ' cle 
relative  to  the  earth  fi:    •        The  vehicle      on  includes  both 
translatioi         bion.  POSITN  also  stores  all  the  traj     ;  ,      ion 
information  in  the  •  ray  POSIT  for  t]     al  output  sccquou.ee  e 
There  is  no  ]  .'  be  3  output  f:        K. 

A  -  2  Jnput •(.'■■■ 

-.,   '  '  •    and  constants  3 
pro;        read  into  i  by  j     read  stai  1  lists 

the  input  variables  and  the  fc    b:  under  ]  '     :y  are  re 

A  r  ber  run  must  be         '         h  vehicle.  If , 

i    .   ;  i         '   be  before  the  vehicle  ]       bhs  surface,  t] 
punched  output  from  the  program  gj    '   '    t  necessary  to  cc 
trajectory  on  the  i   !  run. 

The  ]    ram  out;  '     *s  a  time  history  of  t!   i  I   '  .tion.  The 
dai  a  is  storec  dm  dng  e-    bion  and 

run  by  the  main  p:  ,    ,  D  '   •       execute   '    IST/IECOND/  is 
printed  at  the  c  '    '        This  '       to  aid  the  user-  in 

l°c     ■' '        of  any  ]  termination.  Pre: 

mination  is  generally  cause'  by  j 

The  dimensions  of  the  st«     array  restricts  t!         bo  500 
I  i  i   inter.  Is,  h<  ..ever,  the  program  si  bhal 

■  '  :'        to  K  i0.  Tfcd  s,  (   pled 


;  trajector;  '      '  .  i        un- 

during  a  trajecto] . 
For  :   .'    bput  see  Appendix  C , 


TABLE   1 
]  r  PROGRAM 

Columns   Format    Symbol    Description 

2-8  This  :        ■  '  1  tal '   :■  :<re  than  one  card  i 

includes  the  variables  described  in  Table  2. 
The  list  begins  with  $INCOND  and  ends  ■ 
$SKD. 

2-7  $FLUID    This  namelist  trill  take  one  c      "'  includes 

the  variables  d  '  >le  2C 

2-8  $VLaITS'   This  namelist  will  take  one  card  and  includes 

the  variables  c       1  in  Table  2„ 

2-8  $V:       Th'      list  will  take  one  card  an     Ludes 

the  variables  desc  '    in  1  2. 

2-8  $CONTRL   This  namelist  will  take  one  card  and  incli 

i       variables  des  '  2t 

2-8  $C<       This  i   ' 11  take  mo:       one 

includes  the  variables  lis  bed  in  Tal 

3. 

1  -  10       F  10.0  XI  Th     ]      '  "  to 

11-20      F  10,0  Yl  the  origin  linates. 

21  -  30  0.0  ZI 


Colu  Fornrit         SymT:  .  '  '■.on 

31-40       F  10.0         1    ASS  Tb-    rate  at  ■  particular  ball? 

weight  is  r«  '  .         ds  per  second) , 

I4)    ..  15  IK/  The  number  of  tine  steps  d  ]  ich  a 

toved.     Ti 
is  one  card  oi    I 
from  ■. 
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A  -  3    ]  ' 

'i  j  i     f  olli     Lng  1 

:  2 

LISTS 

Name  List/INCOND/ 


PROGRj  D   UAL 

VARIi  -         JCIATURE       DEFINITION  8    ' 

ISTEP  The  step  numbc 


DT  At  The  time  interval  u     d  d     in;    deba3": 

(seconds) 

T.i  £     t  The  time  elapsed  since  ;    I  <    meed  ( 


i  :       x„  The  x  c    nent  of  the  disi     i  ' 

to  the  inertial  origin  (feet) 


Y£       y„  The  y  coraponenl   E  the  distance  traveled  relative 

to  the  inertial  origin  (fe 

ZE       z  The  z  component  of  the  dis  '.    tr,    i  relatj 

XL 

to  the  inertial  origj      ( 
<f>  Th  le  of  roll  (  ) 

!TA  0  Th     .    g]      of  pitc]    (  is) 

i;'1  IP"  The  angle  of  ya      ( 


'. 


I  I  'AL 

I  ■:  W  ■    '      DE  '] 

s 

:TH    x  The  vehicle  vsd  it;  in 

coord'   '   (feel  per  seco: 

. — s. 

* 

;TH    0  The  vehicle   .  '    :  .■'  .  ,  i  ":     '   .'  ' ,  5 

coord*      (radial   per  seco] 

BVE3      V  The  vehicle    tor  of  translation  ]   GLocity  in 

the  Body  axis  system  (f :   ;  r  sccc. 

BROT     '\  The  vehicle  vector  o::  in  the  Body 

axis  system  (rad: 

W  The  vehicl(    i-ght  (pounds) 

XG       yr  The  x  components  of  the  CG  vector  (feet, 

IG       y  The  y  component  of  t     vector  (fee.': 

ZG       2  The  z   component  of  the  CG  vector  (feet) 

It 

B  B  The  vehicle  bu  ;   '       (]    's) 

XB  xfi  The  x  component  of  t      .  ctor  (feet) 

^3  y^  The  y  c;      '    I   C  vecto 

,/]'  The  7   component  of  1  \        t) 

IXX  I  The  '    ineri  Is  (slugs) 

- 

1  ■  t  of  inert:! 


_i  ',  ■  ■       Ls  (slugs) 
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%  ■  '  Lght  (pounds  ] 

f.  ';  density  ( slugs     < 

:       Z  '  ■   l  ■  ■  '       ; 
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Bpj  .  ]  ■  ■  / 

: '    (; 

Ih©  number  of  ]  Ls  ro- 

les        ' 
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!    v 

IS     iii 

■     i  ' 
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List/Ci 

PROG]  USU. 

\       .  DE    ■ 

:  Xw  The  X  componeni  bor  forces  (pounds) 

1'  The  X  co  .  (         ids) 

SFF 

:  Z  The  Z  component  o:  r  i 

LFF 

I  K_  The  K  c         nent  of    :"  bor  fo:  ids) 

3  The  M  co  b  of  (  ds) 

iiF 

NEFP  The  N  con  t  of  t!  '  .      nds) 

Jir  F 

]  .         List/CO 

X  X  Longitudinal  forces  coeff:  a 


XA  XA.  .  Lor  .     "'    I  <  ' 

novii''  ,        '    rn 


X.  .  Lateral  fo:  ce  coefficients  fc         vehicle 
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YA       XA. .         Lateral  i     •   :    "  nts  f<       cle  movin? 

astern 


Z        Z.  .  No:   ]    ce  coefficients 

ah( 


ZA       2. ' ,.  .. 

as  I 
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movj 


i  '  j    ..  .  Pitch  ]  ',  fore  a  vel  icle 


1"! 


bs  for  a  vehicle  i  lov- 
ins  aste: 


Lpts  i,  j  refer  to  the  possible 
body  vclociti     3  accel        as  d  ble  3« 
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